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Abstract 

By  performing  Quadrature  Mirror  Filter  Bank  (QMFB)  processing  with  a  given 
signal  it  is  possible  to  obtain  Frequency-Time  (F-T)  outputs  that  represent  signal  features 
such  as  bandwidth  (W),  center  frequency  (f  ),  signal  duration  (Tv),  modulation  type  (AM, 
FM,  BPSK,  QAM,  etc),  frequency  content  and  time  allocation.  Because  of  its  unique 
structure,  two  widely  used  signals  based  on  Orthogonal  Frequency  Division  Multiplexing 
(OFDM)  were  chosen  as  signals  of  interest  for  demonstration.  The  general 
implementation  of  the  QMFB  process  is  described  along  with  the  basic  structure  of 
OFDM  signals  related  to  the  physical  layer  perspective  of  802.11a  Wi-Fi  and  802. 16e 
WiMAX  frame  structures  are  described. 

The  adopted  methodology  is  aimed  at  exploiting  signal  of  interest  features 
accounting  for  the  effects  of  signal  resampling  and  zero-padding.  Computed  simulation 
results  are  obtained  after  applying  the  defined  methodology  to  each  signal  of  interest. 
Initial  time  domain  and  frequency  domain  responses  are  presented  for  each  input  signal 
along  with  the  initial  and  computed  resampled  parameters  for  each  case.  Results  for 
selected  QMFB  outputs  are  presented  using  2D  F-T  QMFB  plots  and  ID  average 
frequency  and  average  time  plots.  These  plots  enable  qualitative  visual  assessment  such 
as  may  be  used  by  a  human  operator.  The  ID  responses  are  computed  for  the  input  signal 
and  output  QMFB  responses  and  compared  using  overlay  plots  for  single  burst  and 
multiple  integrated  burst  inputs.  Resultant  time  (At)  and  frequency  (A/)  resolutions  were 
consistent  and  validate  the  usefulness  of  QMFB  processing. 
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OFDM-BASED  SIGNAL  EXPLOT ATION  USING 
QUADRATURE  MIRROR  FILTER  BANK  (QMFB)  PROCESSING 

CHAPTER  1. 

Introduction 

This  chapter  presents  the  research  motivation,  research  objectives  and  research 
organization.  The  research  motivation  is  divided  in  two  subsections  aimed  to  provide  the 
operational  motivation  (Section  1.1.1)  and  the  technical  motivation  (Section  1.1.2)  for  the 
research  effort.  Research  objectives  are  defined  in  Section  1.2.3  with  a  goal  of  1)  finding 
empirical  results  aimed  at  satisfying  the  research  motivation,  and  2)  finding  a  graphical 
representation  that  is  useful  for  highlighting  discriminating  signal  features  from  an 
operators’  perspective.  Finally,  the  research  organization  is  presented  in  Section  1.2.4. 

1 . 1  Research  Motivation 

1.1.1  Operational  Motivation 

Previous  related  work  with  Quadrature  Mirror  Filter  Bank  (QMFB)  processing  [1, 
3,  4]  has  demonstrated  some  practical  capability  for  exploiting  a  given  signal  by  using 
resultant  frequency-time  (F-T)  plots  to  highlight  some  signal’s  distinctive  characteristics. 
Other  signal  exploitation  procedures  using  passive  methods  are  given  in  [8],  which  is 
based  on  performing  wavelet-based  radio  frequency  (RF)  fingerprinting,  and  [0,  11] 
which  is  based  on  performing  Gabor-Based  RF  Distinct  Native  Attribute  (DNA) 
fingerprinting. 

Results  in  previous  works  [1,  3,  4]  are  MATLAB  based  ®  simulated  signals  [1] 
using  the  QMFB  process  to  evaluate  frequency  modulated  CW  (FMCW)  and  binary 
phase  shift  keying  (BPSK)  signals.  Additional  work  has  been  done  with  laboratory  based 
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signals  [4]  using  the  QMFB  process  to  evaluate  noise  technology  radar  (NTR)  signals. 

The  802.11a  Wi-Fi  and  802. 16e  WiMAX  RF  communication  protocols  [9,  12] 
are  continuously  evolving  to  provide  greater  reliability  using  orthogonal  frequency 
division  multiplexing  (OFDM)  techniques  to  better  exploit  available  communication 
channel  resources.  Besides  the  synchronization  and  other  preset  component  features, 
these  signals  present  a  wideband  noise-like  behavior.  This  is  a  result  of  random  symbol 
assignment  making  detection  using  burst  integration  methods  more  difficult.  This 
provides  the  motivation  to  check  if  QMFB  processing  is  applicable  to  experimental 
OFDM-based  signals,  like  802.11a  Wi-Fi  and  802. 16e  WiMAX,  with  an  aim  to 
detennining  if  useful  infonnation  exists  in  QMFB  output  response  that  relates  to  specific 
signal  features. 

1.1.2  Technical  Motivation 

The  technical  motivation  is  aimed  at  presenting  a  new  approach  for  signal 
exploitation  of  OFDM-based  signals.  These  signals  were  chosen  because  of  their  noise¬ 
like  behavior,  and  specifically  the  802.11a  Wi-Fi  and  802. 16e  WiMAX  signals  because 
of  the  randomness  within  signal  generation. 

The  QMFB  process  is  the  baseline  for  conducting  this  research  and  when 
implemented  according  to  its  definition  [1,  3,  4]  has  demonstrated  a  consistent  approach 
to  estimating  signal  parameters  such  as  bandwidth  (W),  center  frequency  (/(),  signal 
duration  (IQ,  modulation  type  (AM,  FM,  BPSK,  QAM,  etc),  frequency  content  and  time 
allocation. 

In  this  research  the  signal  structure  for  each  protocol,  specifically  the  signals’ 


2 


physical  layer,  contains  details  related  to  transmitted  signal  characteristics,  and  those 
details  define  what  the  received  signal  structure  should  be.  By  exploiting  the  QMFB 
process  using  an  OFDM-based  signal  input,  the  goal  is  to  realize  the  extent  that  signal 
characteristics  can  be  reliably  extracted  and  how  well  the  extracted  features  match  the 
defined  structure  [9,  12]. 

1.2  Research  Objectives 

The  main  objective  for  this  research  is  to  perform  a  qualitative  visual  assessment 
of  OFDM-based  signal  responses  using  passive  QMFB  detection.  This  objective  was 
divided  in  two  parts,  including:  1)  reducing  the  required  computation  time  used  in 
previous  works  [1,  3,  4]  to  perform  QMFB  processing,  and  2)  improving  output  results 
related  to  frequency  and  time  resolution  for  a  given  signal  of  interest.  The  goal  is  to 
achieve  frequency  and  time  resolution  that  pennits  reliable  visual  assessment  to  enable 
exploitation  of  a  given  signal’s  characteristics  from  an  operator’s  perspective  within  an 
reduced  computation  time.  This  empirical  approach  would  pennit  exploitation  of  signal 
characteristics  such  as  bandwidth  (W),  center  frequency  (f),  signal  duration  (Tv), 
frequency  content  and  time  allocation  as  presented  in  2D  QMFB  F-T  plots. 

1 .3  Research  Organization 

The  document  includes  general  descriptions  and  infonnation  for  specific  cases 
that  were  used  to  compute  the  results  according  to  the  research  objectives  defined  in 
previous  section. 

Chapter  2  presents  the  necessary  technical  background  used  as  a  baseline  during 
this  research  effort.  Description  about  QMFB  process  is  presented  along  with  an  OFDM 
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overview  and  some  basic  characteristics  of  the  signals  of  interest  with  emphasis  on 
physical  layer  characteristics  [9,  12]. 

Chapter  3  presents  the  adopted  methodology  aimed  to  compute  the  necessary 
results  to  achieve  the  defined  objectives  according  to  the  technical  background  defined  in 
previous  chapter.  A  process  overview  is  presented  and  its  decomposition  is  described. 
The  post-collection  process  is  defined  first,  followed  by  verification  and  validation  of  the 
QMFB  process.  The  input  signal  parameters  are  then  verified  according  to  the  standards 
defined  for  each  case  [9,  12].  The  QMFB  process  is  described,  including  the  effects  of 
signal  resampling  and  zero-padding  used  for  given  input  signals.  Measurable  process 
outputs  are  defined,  signal  of  interest  parameters  are  provided,  and  the  graphical 
presentation  fonnat  for  QMFB  results  is  introduced. 

Chapter  4  presents  computed  simulation  results  according  to  the  given 
methodology  and  comparisons  to  initial  parameters  for  each  signal  of  interest  are 
presented.  Initial  conclusions  for  each  case  are  also  presented.  The  results  are  presented 
in  2D  F-T  plots  along  with  corresponding  ID  plots  for  time  domain  and  frequency 
domain  responses.  Overlay  plots  are  used  for  initial  input  signal,  output  signal  and  burst 
integrated  output  signal  comparisons. 

Chapter  5  presents  the  research  summary  and  conclusion.  The  motivation, 
methodology  and  computed  versus  expected  results  are  discussed.  Recommendations  for 
future  work  are  given.  Finally,  appendices  are  provided  for  each  signal  of  interest  along 
with  the  implemented  MATLAB®  code. 
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CHAPTER  2. 


Background 

2.1  Overview 

This  chapter  provides  the  technical  background  established  for  this  research 
effort.  This  chapter  is  divided  in  two  main  sections  aimed  at  defining  a  performance 
baseline  and  providing  basic  concepts  related  to  the  topics  and  experimental  techniques 
exploited  during  this  research.  The  Quadrature  Mirror  Filter  Bank  (QMFB)  process  is 
described  in  Section  2.2,  which  provides  general  implementation  parameters  and 
assumptions.  Fundamentals  of  Orthogonal  Frequency  Division  Multiplexing  (OFDM)  are 
introduced  in  Section  2.3  which  provides  the  general  signal  structure  as  well  as  key  signal 
features.  There  were  two  OFDM-based  signals  of  interest  for  this  research.  The  physical 
layer  characteristics  of  each  are  provided  in  Section  2.4  for  the  802.11a  Wi-Fi  signal  and 
Section  2.5  for  the  802. 16e  WiMAX  signal. 

2.2  Quadrature  Mirror  Filter  Bank  (QMFB)  Processing 

The  QMFB  is  an  orthogonal  waveform  decomposition  technique  based  on 
wavelet  filter  theory.  Each  layer  output  provides  input  signal  frequency  and  time 
characteristics  aimed  to  estimate  and  exploit  various  signal’s  features.  Common  features 
of  interest  include  signal  modulation  type,  bandwidth,  frequency  component  distribution, 
signal  duration  as  well  as  time  and  frequency  allocation.  As  detailed  in  Figure  2.1,  the 
process  adopted  here  was  introduced  by  Pace  [1]  and  subsequently  exploited  in  additional 
related  works  of  Jarpa  and  Atienza  [3,  4],  Jarpa’s  research  in  [3]  was  based  on  verifying 
QMFB  response  given  structured  signal  inputs,  including  a  single  tone,  multiple  tone, 
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frequency  modulated  continuous  wave  (FMCW),  and  frequency/phase  shift  keyed 
(FSK/PSK)  waveforms.  Atienza’s  research  in  [3,  4]  analyzed  the  QMFB  response  given 
noise-like  input  signals  collected  from  AFIT’s  Noise  Technology  Radar  (NTR)  Noise 
Network  (NoNET). 

For  a  given  input  signal,  the  QMFB  process  performs  waveform  decomposition  by 
computing  continuous  signal  decimation  in  high  and  lower  frequency  components  using 
the  signal  structure  detailed  in  Figure  2.1.  The  input  signal  is  decomposed  into  in-phase 
(I)  and  quadrature  (Q)  components  using  modified  sine-shaped  filters  [1]. 


Layer  1  Layer  2  Layer  3 


Tjme  Time  Time  Time 

Figure  2.1  QMFB  process  overview  [1] 
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The  output  of  each  filter  “Layer”  corresponds  to  the  amplitude  (real  or  complex) 
or  magnitude  as  a  function  Frequency-Time  (F-T)  parameters  according  to  the  input 
signal.  After  each  filtering  stage,  the  signal  is  decimated  so  further  layers  can  be 
computed.  Because  of  the  decimation  process,  and  the  filter  input  signal  being  a  function 
of  the  previous  layer,  the  initial  time  window  (duration)  is  increased  by  a  power  of  two 
with  every  computation.  On  the  other  hand,  the  initial  frequency  window  (bandwidth) 
decreases  according  to  the  same  power  of  two.  This  creates  a  tradeoff  between  the 
different  layers  which  could  produce  useful  data  for  a  given  layer.  Given  the  post¬ 
filtering  decimation,  the  total  numbers  of  available  layers  (N)  is  a  function  of  the  initial 
number  of  data  samples  (Ns)  and  given  by: 

N=log2(N,)  (2.1) 

For  a  given  N,  the  frequency  window  extent  (Af)  and  time  window  duration  (At) 
for  a  given  layer  (N£)  can  be  computed  as  follows: 


A/ 


fs 


2(2N~l  -1) 


(2.2) 


At  = 


fs(2N  ‘  ~  1)  fsNT 


(2.3) 


Where  fs  is  the  sampling  frequency,  Nf  is  a  given  frequency  layer  number,  and  Nr 
is  a  given  time  layer  number.  According  to  (2.2)  and  (2.3)  it  can  be  seen  that  lower 
numbered  layers  present  shorter  time  windows  compared  to  higher  layers  and  there  is  an 
inverse  relationship  with  the  frequency  window. 

The  decimation  is  computed  according  to  [1],  which  was  proven  to  work  in 
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previous  works  [3,  4]  and  it  consist  in  the  implementation  of  a  “modified  sine”  finite 
impulse  response  (FIR)  given  by: 


,  r  7  S  .  ( n+0.5 
h  n /=.  —  sine  - 

Vi  {  c 


N  (N-2) 

win], - <  n  <  - - 

2  2 


Where  S  is  the  scaling  variable,  C  is  the  compression  variable,  N  is  the  number  of 
coefficients,  and  w[n]  is  a  Hamming  window  weighting. 

These  particular  filters  have  a  flat  bandpass  response  and  pass  the  maximum 
amount  of  desired  signal  energy  at  each  layer.  According  to  [1,  3]  with  an  aim  of  getting 
“nearly  orthogonal  filters  with  cross-correlation  of  less  than  0.001”  using  N=  512,  the 
constant  compression  and  scaling  values  in  (2.4)  are  C=l. 99375872338059  and 
S=l. 00618488680080.  The  Hamming  window  is  use  to  suppress  the  effects  of  Gibb’s 
phenomena  resulting  from  sequence  truncation  [1,3]. 

To  avoid  data  sample  loss,  and  because  the  total  number  of  available  layers  is 
given  by  (2.1),  the  following  assumption  has  been  made  to  compute,  as  initial 
approximation,  the  total  number  of  available  layers  for  a  given  input  signal: 

N=c  eil[log2  (Ns )]  (2.5) 


Where  N  is  the  total  number  of  available  layer,  Ns  is  the  number  of  samples,  and 
the  ceil[  ]  operation  ensures  zero-padding  to  the  next  power  of  two. 

2.3  Orthogonal  Frequency  Division  Multiplexing  (OFDM)  signals 

The  sub-carrier  frequencies  in  OFDM  are  chosen  to  be  mutually  orthogonal  and 
inter-carrier  guard  bands  are  not  required  as  in  basic  modulation  process.  This  simplifies 
the  transmitter  and  receiver  designs;  unlike  conventional  FDM,  which  requires  a  separate 
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filter  for  each  sub-channel. 


Considering  the  frequency  domain,  an  OFDM  symbol  is  made  up  multiple 
subcarriers,  the  number  of  which  determines  the  Fast  Fourier  Transform  (FFT)  size  used. 
As  illustrated  in  Figure  2.2  [12]  there  are  three  distinct  types  of  subcarriers  used.  The 
type  of  subcarrier  and  purpose  are  as  follows: 

•  Data  Subcarriers:  Data  Transmission 

•  Pilot  Subcarriers:  Signal  Estimation 

•  Null  Carriers:  No  transmission;  guard  bands  and  DC  subcarrier. 

Data  Subcarricrs  DC  subcarrier  Pilot  subcarriers 


Figure  2.2  OFDM  frequency  structure  [12] 

Considering  time  domain,  an  Inverse  Fourier  Transform  (IFT)  creates  the  OFDM 
wavefonn  where  the  signal  time  duration  ( Ts)  is  the  result  of  the  initial  Guard  Time  ( 7V) 
plus  the  useful  symbol  time  ( Tb ).  A  Tg  copy  of  the  last  useful  symbol  is  added  and  used 
to  correct  for  multipath.  Therefore,  the  basic  OFDM  time  structure  is  given  by  Figure 
2.3. 
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Ts 


Figure  2.3  OFDM  time  structure  [12] 


The  resultant  transmitted  signal  is  given  by 

/ 


s(t)=  Re 


Nused 

2 


eW  ^  ckd2nkA^ 


ifc-Tg) 


\ 


(2.6) 


kyO 


Where  fc  is  the  center  carrier  frequency,  Nuseci  is  the  number  of  used  subcarriers,  cy 
is  a  complex  modulation  number  specifying  a  point  in  the  QAM  signaling  constellation, 
A f  is  the  subcarrier  spacing,  t  is  the  elapsed  time  since  the  beginning  of  the  OFDM 
symbol,  and  Tg  =  GxTh  with  defined  G  e  [1/4,  1/8,  1/16,  1/32]  [12]. 

According  to  the  OFDM  frequency  and  time  definitions  given  above,  an  OFDM- 
based  signal  is  generated  as  given  by  Figure  2.4 


Figure  2.4  OFDM-Based  Signal  Generation  Process 
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2.3.1  OFDM-Based  802. 1  la  Wi-Fi  Signal 

This  OFDM-based  signal  is  widely  used  in  the  implementation  of  wireless  local 
area  networks  (WLAN).  According  to  [9]  and  related  to  this  research  effort,  the  OFDM 
802.11a  signal  covers  the  following  frequencies  5.15-5.25  GHz,  5.25-5.35  GHz,  and 
5.725-5.825  GHz.  This  provides  a  wireless  LAN  with  data  payload  communication  bit 
rates  of  Rb  e  [6,  9,  12,  18,  24,  36,  48,  54]  Mbit/s.  The  system  uses  Ns  =  52  subcarriers 
that  are  modulated  using  binary  phase  shift  keying  (BPSK),  quadrature  phase  shift  keying 
(QPSK)  or  16-ary,  64-ary  quadrature  amplitude  modulation  (16-QAM  or  64-QAM).  The 
defined  802.1  la  Wi-Fi  signal  frame  fonnat  [9]  according  to  the  Physical  Layer  Converge 
Procedure  (PLCP)  is  defined  as  the  PLCP  Protocol  Data  Unit  (PPDU)  and  includes  the 
OFDM  PLCP  preamble,  OFDM  PLCP  header  and  the  PLCP  service  data  unit  (PSDU)  as 
shown  in  Figure  2.5. 


Rate 

Reserved 

LENGTH 

Parity 

Tail 

SERVICE 

PSDU 

Tail 

4  bits 

1  bit 

1 2  bits 

1  bit 

6  bits 

1 6  bits 

6  bits 

(BPSK,  r='A) 


(RATE  is  indicated  in  signal) 


PLCP  Preamble  SIGNAL  DATA 

12  Symbols  One  OFDM  Symbol  Variable  Number  of  OFDM  Symbols 


Figure  2.5  PPDU  frame  fonnat  [9] 

The  PLCP  preamble  field  is  present  for  the  receiver  to  acquire  an  incoming 
OFDM  signal  and  synchronize  the  demodulator.  The  preamble  consists  of  12  total 
symbols,  including:  1)  ten  short  symbols  for  establishing  automatic  gain  control  (AGC), 
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coarse  carrier  frequency  estimation,  and  2)  two  long  symbols  for  fine  frequency 
acquisition  in  the  receiver.  The  PLCP  preamble  structure  is  shown  in  Figure  2.6. 


8  +  8  =  16  fis- 


◄ - 10  x  0.8  =  8  |is - ► 

- 2  x  0.8  +  2  x  3.2  =  8  gs - 

•0.8+3. 2=4  (j,s- 

■0.8+3. 2=4  ns- 

■0.8+3. 2=4  ns- 

ti 

t2 

t3 

u 

t5 

t6 

t7 

B 

t9 

tio 

GI*2  |  T,  J  T2 

GI 

Data  1 

Signal  Detect,  AGC,  Coarse  Freq.  Channel  and  Fine  Frequency  Rate  Service  +  Data 

Diversity  selection  Offset  Estimation  Offset  Estimation  Length 

Timing 
Synchronize 

Figure  2.6  802.1  la  Wi-Fi  Signal  preamble  structure  [9] 

2.3.2  OFDM-Based  802. 16e  WiMAX  Signal 

This  protocol  is  aimed  to  extend  wireless  range  of  previous  WLAN  protocols  and 
provide  broadband  connectivity  for  data  and  telecommunications.  According  to  [12]  and 
related  to  this  research  effort,  the  OFDM  802. 16e  signals  cover  the  frequency  bands 
below  11  GHz.  In  this  case,  each  data  frame  is  divided  in  two  subframes,  including: 
1)  the  down  link  (DL)  subframe  aimed  to  transmitting  data  and  control  messages  to 
specific  subscriber  station  (SS)  and  2)  the  up  link  (UL)  subframe  that  is  used  by  the 
subscriber  to  transmit  to  the  Base  Station  (BS).  Each  subframe  can  be  modulated  using 
BPSK,  QPSK,  16-QAM  or  64-QAM.  The  defined  802. 16e  WiMAX  Signal  frame  can  be 
fonnatted  using  either  time  division  duplexing  (TDD)  or  frequency  division  duplexing 
(FDD)  techniques.  As  shown  in  Figure  2.7  for  TDD,  the  generic  802. 16e  WiMAX  Signal 
frame  time  duration  is  obtained  by  adding  each  subframe  (DL  and  UL)  per  SS. 
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Figure  2.7  OFDM  frame  structure  with  TDD  [12] 

As  shown  in  Figure  2.8  for  FDD,  the  generic  802. 16e  WiMAX  Signal  frame  time 
duration  is  constant.  The  subframes  (DL  and  UL)  for  all  SSs  use  the  same  time  frame  but 
with  different  frequencies. 

Time 


Figure  2.8  OFDM  frame  structure  with  FDD  [12] 
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2.4  Summary 

Technical  background  for  the  research  effort  has  been  presented,  to  include  a 
discussion  of  the  QMBF  process  and  characteristics.  The  total  number  of  samples,  sample 
frequency  and  zero  padding  were  described  as  key  parameters  aimed  to  achieve  Af  and  At 
which  allow  reliable  qualitative  visual  assessment  according  to  the  available  generated 
layers.  Fundamentals  of  OFDM  were  introduced  and  two  specific  OFDM-based  signals 
of  interest  described.  This  included  a  discussion  of  relevant  physical  layer  characteristics 
for  802. 11a  Wi-Fi  and  802. 16e  WiMAX  signals. 
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CHAPTER  3. 


Methodology 

3.1  Introduction 

This  chapter  discusses  the  adopted  methodology  aimed  at  perfonning  the 
necessary  data  and  simulation  management  to  satisfy  defined  objectives  in  this  the 
research  effort.  Section  3.2  provides  the  process  overview  which  shows  the  main  flow 
diagram  that  was  used  as  the  baseline  process  for  the  research  effort.  Section  3.3 
describes  the  verification  and  validation  (V&V)  signals  used  to  ensure  the  QMFB  process 
was  implemented  correctly.  This  is  followed  by  Section  3.4  which  describes  the  OFDM- 
based  signals  considered  for  demonstrating  the  exploitation  capability  of  the  QMFB 
process.  Implementation  of  the  QMFB  process  is  detailed  in  Section  3.5,  to  include  the 
effects  of  signal  resampling  and  zero-padding  for  a  given  input  signal.  Finally,  the 
chapter  concludes  with  Section  3.6,  results  presentation  fonnat,  which  shows  how  results 
are  presented  for  each  case  considered. 

3.2  Process  Overview 

The  flow  diagram  in  Figure  3.1  was  developed  to  set  the  sequence  of  steps  aimed 
at  achieving  the  objectives  described  in  Chapter  1  and  based  on  background  information 
given  in  Chapter  2.  The  goal  was  to  provide  measurable  results  at  every  different  stage 
of  the  modeled  problem.  For  each  input,  signal  characteristics  were  first  verified  using 
both  time  domain  and  frequency  domain  power  spectral  density  (PSD)  responses.  Once 
verified,  the  signal  was  input  to  the  QMFB  process  and  the  resulting  layer  outputs  were 
used  to  exploit  features.  The  exploitation  assessment  included  two  steps:  1)  Seeing  how 
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well  resultant  QMFB  features  match  the  1-D  time  and  PSD  characteristics  of  the  input 
signal,  and  then  2)  Using  qualitative  visual  assessment  to  see  if  the  2-D  QMFB  F-T 
outputs  provided  an  additional  insight  on  features  not  evident  in  1-D  responses. 


YES 


Figure  3 . 1  Overview  of  Research  Methodology 

3.3  Verification  and  Validation  (V&V)  signals 

To  ensure  the  QMFB  process  was  implemented  correctly,  QMFB  output 
responses  were  looked  at  using  two  specific  input  signals  for  V&V.  The  two  analytic 
V&V  signals  included  a  continuous  LFM-modulated  signal  and  a  discrete  multi-tone 
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signal  used  for  V&V  in  previous  work  [1,  3,  4].  The  signal’s  characteristics  presented  on 
each  V&V  signals  are  aimed  to  realize  differences  on  the  QMFB  response  due  to 
continuous  modulation,  discrete  modulation,  and  single  versus  multiple  frequencies. 
Signal  generation  parameters  for  each  of  the  V&V  signal  are  set  to  establish  identical 
QMFB  processing  parameters  as  provided  in  Table  3.1 


Table  3.1  QMFB  Parameters  for  V&V  Signal  Processing 


Bandwidth 

Samp  Freq 

Duration 

Number  of 

(KHz) 

fs  (Hz) 

(mScc) 

Samples 

40 

100 

1.6 

2428 

3.3.1  Continuous  Linear  FM  (LFM)  signal 

The  LFM  V&V  signal  was  used  to  assess  the  QMFB  response  to  a  continuous 
signal  input  having  a  linear  frequency  and  amplitude  change  during  the  signal  duration 
[2],  The  signal  time  and  PSD  response  is  shown  in  Figure  3.2  were  generated  using  the 
analytic  expression  in  (3.1)  with  fL  =  45  KHz ,  /)/=95KHz,  and 

fA  =  (Fl  -  F„)  =  -50  KHz. 


s  |  (/)  =  A  /  .s//j  (27r///t) 

S2  (0  =  A2(t)  sin(2nfhf) 
s3  (t)  =A3sin{2nfLt) 


0  <  t  <  0.6  m 
0.6  m  <  t  <  1.0  m 
1.0  m<  t  <1.6  m 


J 

s/.FA/(0  '  Sj(t) 

i=l 


(3.1) 
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Figure  3.2  Analytic  LFM  V&V  signal  time  domain  and  PSD  responses. 

3.3.2  Discrete  Multi-Tone  (DMT)  V&V  Signal 

The  DMT  V&V  signal  was  generated  consistent  with  Szmajda’s  V&V  signal  in 
[7]  and  was  chosen  to  assess  the  QMFB  response  to  a  discrete  modulated  signal  having 
both  single  and  multiple  frequency  components  across  time.  The  signal  time  and  PSD 
responses  are  shown  in  Figure  3.3  and  were  generated  using  the  analytic  expressions  in 


(3.2)  and  tone  parameters  provided  in  Table  3.2 


Si  (0  = 

=  A 1  sin(27ifjt) 

0 

A 

t 

< 

13 

82(0  = 

=  A2  sin{2nf^t) 

13  < 

t 

< 

25 

83(0  = 

=  0 

25  < 

t 

< 

37 

84(0  = 

=  Sj(t) 

37  < 

t 

< 

50 

85(0  = 

=  A5sin{2irf.i) 

37  < 

t 

< 

50 

SfiCO  = 

=  A6sin{2uf6t ) 

37  < 

t 

< 

50 

u 

sdmt(.o= 'y'  Si® 

i=l 
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Table  3.2  Discrete  Multi-Tone  Szmajda  Signal  Generation  Parameters 


Tone 

Si(t) 

*2(0 

s4(t) 

*5(0 

Sd(0 

Amplitude 

230V2 

2 Aj 

A, 

Aj 

Aj 

Frequency  (Hz) 

5 

5 

5 

10 

40 

Ll_ 

-  -20 


-30 - - ■ - 1 - 1 - 1 - 1 - 1 -  - 

0  5  10  15  20  25  30  35  40  45 

Frequency  (Hz) 

Figure  3.3  Analytic  DMT  Szmajda  V&V  signal  time  domain  and  frequency  responses. 

3.4  OFDM-Based  Signals 

The  QMFB  response  to  the  V&V  signals  described  in  Section  3.2  and  Section  3.3 
provide  a  baseline  for  assessing  QMFB  exploitation  potential  using  two  OFDM-based 
signals.  Each  signal  was  decomposed  according  to  its  specific  characteristics  in  order  to 
find,  evaluate  and  exploit  different  responses.  The  two  signals  considered  include 
experimentally  collected  802.11a  Wi-Fi  [8]  and  802. 16e  WiMAX  [0,  1 1]  signals. 
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3.4.1  802.1  la  Wi-Fi  signal 


This  OFDM  signal  was  chosen  because  of  its  well-known  and  structured  physical 
layer  response  as  described  in  Chapter  2.  In  this  case,  the  signal  used  corresponds  to 
experimental  collections  taken  to  support  results  in  [8].  The  signal  was  analyzed  in  three 
stages  using  1)  an  isolated  preamble  response,  2)  a  single  burst  response  and  3)  an 
integrated  burst  response. 

During  the  first  analysis  stage  the  802.11a  preamble  response  was  isolated  and 
only  the  first  10  short  symbols  corresponding  to  the  first  half  of  the  preamble  was 
considered.  During  the  second  analysis  stage,  the  response  of  a  single  802.1  la  burst  was 
set  as  input  to  the  QMFB  process.  Finally,  burst  integration  was  computed  considering 
the  mean  amplitude  response  for  Nb  =  500  burst  collections.  The  802.11a  Wi-Fi  signal 
parameters  for  each  analysis  case  are  presented  in  Table  3.3 


Table  3.3  802. 11a  Wi-Fi  Signal  Parameters 


Input 

Bandwidth 

(MHz) 

Samp  Freq 
fs  (MHz) 

Duration 

(pSec) 

Number  of 
Samples 

Preamble 

8.0 

23.75 

135 

3200 

Burst 

9.0 

23.75 

124 

2945 

3.4.2  802. 16eWiMAX  signal 

This  signal  was  chosen  to  assess  the  QMFB  response  using  a  more  complex 
structured  OFDM-based  signal  input.  In  this  case,  the  input  signal  corresponds  to 
experimental  collections  made  in  support  of  work  in  [0,  11].  The  process  is  aimed  to 
exploit  the  physical  layer  parameters  and  structure  described  in  Chapter  2  according  to 
[12].  The  experimental  802. 16e  WiMAX  signal  was  analyzed  in  two  stages.  First,  the 
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range-only  mode  response  was  considered  as  input  for  two  cases,  including  single  burst 
and  burst  integration  using  Nb~  1400  bursts.  Secondly,  the  data-only  mode  response  was 
considered  as  input  for  two  cases,  including  single  burst  and  burst  integration  using 
Nil  ~  640  bursts.  The  802. 16e  WiMAX  signal  parameters  for  each  case  described  above 
are  presented  in  Table  3.4 


Table  3.4  802. 16e  WiMAX  Signal  Parameters 


Mode 

Bandwidth 

(MHz) 

Samp  Freq 
fs  (MHz) 

Duration 

(mScc) 

Number  of 
Samples 

Range-Only 

2.0 

11.875 

0.381 

4521 

Data-Only 

2.5 

11.875 

1.61 

19125 

3.5  QMFB  Processing 

Aimed  to  evaluate  experimental  results  according  to  the  available  inputs  described 
in  Section  3.3  and  Section  3.4,  the  QMFB  process  was  first  modeled  in  MATLAB®  based 
on  previous  work  [1,  3,  4]  and  according  to  the  process  described  in  Chapter  2.  The 
length  of  each  input  signal  vector  was  first  checked  and  zero-padded  to  the  closest  integer 
power  of  2  to  avoid  data  loss  during  QMFB  processing  (2.5).  Then,  the  first  simulation 
trial  was  completed.  After  verifying  that  the  QMFB  process  was  working  for  every 
available  input  signal  type,  multiple  trials  were  completed  to  check  differences  and  to 
choose  the  best  QMFB  layer  output  response  according  to  the  specific  input.  According 
to  the  available  data  format  for  a  given  input,  up  to  four  different  cases  were  used  in  the 
modeled  QMFB  process.  The  initial  simulation  parameters  used  to  check  QMFB 
performance  are  shown  in  Table  3.5 
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Table  3.5.  Initial  QMFB  Configuration 


Input 

Real  (R),  Complex(C),  Magnitude 

#  Filter  Coefficients 

512 

Window 

Hamming 

fs 

Initial  Sample  Frequency 

Total  #  QMFB  Layers 

Ceil  { logo  [length(  input)] } 

Given  an  input  vector  to  the  QMFB  process,  the  first  “Layer”  output  is  computed, 
decimated,  stored,  and  passed  to  the  next  layer.  This  iterative  process  repeats  and  ends 
after  the  initial  vector  is  decimated  according  to: 

Nl  =log2  ( Data  Length )  =Total  #  Available  Layers  (3.3) 

Because  data  length  is  a  function  of  both  sampling  frequency  (fs)  and  total  signal 

duration  (Ts),  two  pre-conditioning  steps  were  included  in  the  QMFB  process,  henceforth 
referred  to  as  resampling  and  zero-padding.  Thus,  each  input  signal  was  resampled,  zero- 
padded  and  set  as  input  to  the  QMFB  process  to  achieve  the  desired  frequency  and  time 
resolution  in  a  given  QMFB  layer  output.  According  to  the  QMFB  frequency-time 
tradeoff  described  in  Chapter  2,  better  time  resolution  is  achieved  in  lower  QMFB  layers, 
but  better  frequency  resolution  is  achieved  in  higher  QMFB  layers.  So  the  challenge  is  to 
find  a  frequency  sampling  rate  ( Q  versus  signal  duration  (Ts)  aimed  to  present  an 
accurate  representation  of  input  signal  characteristics.  Henceforth,  this  is  called  the  “the 
most  representative  layer”  which  corresponds  to  the  layer  which  shows  frequency  and 
signal  time  allocation  that  could  be  useful  to  exploit  input  signal  features  such  as  signal 
bandwidth  (W),  signal  duration  (Ts),  time  resolution  (At)  and  frequency  resolution  (Af). 
Figure  3.4  shows  a  collection  of  QMFB  outputs  where  “the  most  representative  layer” 
corresponds  to  layer  #20  (Q20)  highlighted  in  the  red  dashed  rectangle. 
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Frequency  [Hz]  Frequency  |Hz|  Frequency  [Hz]  Frequency  [Hz] 


Figure  3.4  QMFB  output  collection  sample 
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3.5.1  Signal  Resampling 

By  resampling  the  input  signal,  the  number  of  available  input  data  samples 
increases  without  changing  the  total  signal  duration  Ts.  This  also  increases  the  effective 
sampling  frequency  f  and  because  of  the  F-T  trade-off  in  QMFB  processing,  it  is 
expected  that  “the  most  representative  layer”  should  be  located  in  upper  layers,  which 
increases  required  computing  time  and  time  resolution  At  for  a  given  input  signal. 

3.5.2  Zero  Padding 

By  zero-padding  the  input  signal,  the  total  effective  time  duration  Ts  is  increased, 
so  better  frequency  resolution  is  achieved  for  a  given  signal  and  it  is  expected  that  the 
most  representative  layer  is  now  located  in  lower  layers.  But  lower  layers  present  poorer 
signal  frequency  resolution  (increased  Af),  therefore  by  using  lower  layer  analysis  results 
could  be  an  inaccurate  signal  representation. 

3.5.3  Measurable  outputs 

Aimed  to  extract  useful  data  and  to  reduce  computation  time  of  the  QMFB 
process,  the  QMFB  process  was  divided  in  four  steps,  named  as:  1)  layer  generation 
minimization,  2)  Qn  matrix  formatting,  3)  time  analysis  and  4)  frequency  analysis.  By 
perfonning  Layer  Generation  Minimization  it  was  possible  to  reduce  the  total  number  of 
generated  layers,  because  higher  layers  does  not  improve  qualitative  visual  analysis  and 
can  be  disregarded.  The  Qn  matrix  formatting  is  a  matrix  dimensional  reduction  aimed  to 
disregard  irrelevant  analysis  data  such  as  induced  by  zero-padding  and  high  frequency 
data  falling  outside  the  processed  bandwidth  of  the  input  signal.  In  this  case,  the  first 
useful  plots  can  be  computed  considering  three  variables,  including  1)  desired  layer 
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number,  2)  upper  frequency  limit,  and  3)  upper  time  limit.  The  resultant  output  Qn.mat 
file  was  created  (were  n  denotes  a  given  layer  number)  considering  an  amplitude  (real  or 
complex)  matrix,  a  magnitude  matrix,  At  and  Af  values  for  a  given  layer,  time  vector, 
frequency  vector  and  fs.  The  time  analysis  and  frequency  analysis  steps  permit 
comparison  of  input  data  and  QMFB  responses  in  order  to  realize  process  accuracy  and 
losses  due  to  the  signal  processing.  Figure  3.5  shows  all  the  available  measurable  outputs 
and  their  relation  to  the  QMFB  computing  process. 
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Figure  3.5  Measurable  process  outputs 


3.5.4  Presentation  of  QMFB  Layer  Outputs 

Aimed  to  the  objective  of  providing  QMFB  qualitative  visual  assessment,  the 
results  are  presented  in  frequency  versus  time  (F-T)  plots  for  “the  most  representative 
layer”  along  with  the  average  frequency  and  average  time  plot  of  the  input  processed 
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signal.  The  QMFB  layer  output  data  can  be  presented  in  many  formats.  For  clarity  and 
to  enable  consistent  comparison  as  the  input  signal  varies,  all  results  presented  in 
Chapter  4  use  the  format  presented  in  Figure  3.6  which  includes  numbered  responses  of: 

1 .  The  normalized  average  frequency  response  computed  as  a  row-wise  average  for 
a  given  QMFB  layer  data. 

2.  “The  most  representative  layer”,  corresponds  to  a  2D  QMFB  F-T  layer  output 
plotted  using  Mat  lab"  pcolor  function  followed  by  shading  interp.  Appendix  B 
contains  a  complete  set  of  2D  F-T  for  given  signal. 

3.  The  normalized  average  time  response  computed  as  a  column-wise  average  for  a 
given  QMFB  layer  data. 


Figure  3.6.  Representative  Presentation  Layout  for  a  given  QMFB  layer  output  showing 
(1)  1-D  Average  Frequency,  (2)  2-D  QMFB,  and  1-D  Average  Time  responses 

Note:  For  Chapter  4  results,  the  plots  correspond  to  three  different  QMFB  layers 
identified  as  the  most  representative  responses  for  a  given  domain. 
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3.5.5  Single  Burst  vs.  Integrated  Bursts  Overlay  Plots 

Aimed  to  check  QMFB  response  related  to  a  given  signal  input,  average  plots  for 
both  time  (Figure  3.7)  and  PSD  (Figure  3.8)  responses  are  provided  to  assess  process 
accuracy  and  signal  processing  loses.  The  following  two  plots  are  representative  of  the 
presentation  format  used  in  Chapter  4. 


Figure  3.7  Overlay  Time  responses  showing  input  signal  response,  QMFB  output 

for  single  and  integrated  burst 
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Figure  3.8  Overlay  Frequency  responses  showing  input  signal  response,  QMFB 
output  for  single  and  integrated  burst 


3.6  Summary 

The  methodology  presented  in  this  chapter  was  applied  to  each  of  the  described 
input  signals.  To  reduce  computation  time,  the  QMFB  process  was  decomposed  and 
measurable  outputs  at  each  stage  were  defined.  Using  the  designed  process  flow  diagram 
in  Figure  3.5,  different  layers  output  data  are  generated,  saved  and  analyzed  using 
qualitative  visual  assessment  to  characterize  QMFB  potential  for  exploiting  unknown 
signals.  The  effectiveness  of  this  method  is  based  on  signal  parameters  such  as  the 
resampling  vector,  zero-padding  factor,  total  signal  duration,  time  resolution  At  and 
frequency  resolution  Af  for  each  computed  QMFB  layer. 


28 


CHAPTER  4. 


Simulation  Results  and  Analysis 

4.1  Introduction 

This  chapter  presents  MATLAB®  simulation  results  and  data  analysis  based  on 
that  was  obtained  using  the  methodology  discussed  in  Chapter  3.  Baseline  verification 
and  validation  (V&V)  perfonnance  of  the  QMFB  process  is  first  addressed  in  Section  4.2 
using  the  LFM  and  analytic  Szmajda  signals  described  in  Chapter  2.  Then,  OFDM  based 
signals  described  in  Chapter  2  are  input  to  the  QMFB  process  to  determine  if  visually 
discernible  features  are  present  for  estimating  signal  parameters.  The  QMFB  results  are 
presented  in  frequency  versus  time  (F-T)  plots  for  a  given  layer  previously  defined  in 
chapter  3  as  “the  most  representative  layer”,  along  with  individual  average  time  and 
frequency  responses.  Finally,  overlaid  plots  are  presented  to  compare  QMFB  process 
outputs  for  single  burst  and  integrated  burst  responses. 

4.2  Process  V&V  Signals 

Perfonnance  assessment  is  first  performed  with  the  LFM  and  analytic  Szmajda 
signals  described  in  Chapter  2  input  to  the  QMFB  process.  QMFB  perfonnance  is 
characterized  trough  qualitative  visual  assessment  using  joint  2D  F-T  responses  as  well 
as  independent  ID  frequency  and  time  responses. 

4.2. 1  LFM-Based  V&V  Signal 

The  normalized  time  and  frequency  responses  for  the  input  LFM-based  signal  are 
presented  in  Figure  4.1.  The  signal  was  resampled  and  zero-padded  prior  to  QMFB 
processing  according  to  the  values  shown  in  Table  4.1.  The  LFM  signal  time  response 
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shows  that  higher  frequencies  are  located  in  the  first  half  of  the  signal,  after  those  lower 
frequency  responses  are  present.  Related  to  the  PSD  response,  it  can  be  seen  that  the 
WsodB  bandwidth  is  located  within  the  /=  40  KHz  to  /=  100  KHz  range.  Two 
frequencies, /=  45  KHz  and/=  95  KHz,  are  present  with  higher  power  levels  of  P  =  0  dB 
and  P  =  -4  dB,  respectively.  For  the  rest  of  the  frequencies  between  the  two  power  peaks 
it  can  be  seen  that  there  is  an  inverse  relationship  between  the  frequency  and  signal 
power,  with  power  decreasing  from  P  =  -16  dB  to  P  =  -19  dB. 


Figure  4. 1  LFM  V&V  signal  time  domain  and  PSD  responses. 


Table  4.1  LFM  V&V  Signal  Parameters 


Input 

Bandwidth 

(Hz) 

Samp  Freq 
fs  (Hz) 

Duration 

(Sec) 

Number  of 
Samples 

Sample 

Rate 

Zero 

Padding 

Original 

6.0  x  104 

1.5  x  108 

1.62  x  10'3 

2428 

1 

N/A 

Resampled 

6.0  x  104 

1.2  x  1011 

3.50  x  10'3 

1942400 

800 

222 
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After  pre-processing  the  LFM  signal  according  to  parameters  in  Table  4.1,  the 
signal  was  input  to  the  QMFB  process.  Representative  QMFB  results  for  Layer  #14  are 
presented  in  Figure  4.2  with  the  signal’s  linear  frequency  behavior  highlighted  (yellow 
arrows)  and  power  distributed  according  to  color  bar.  The  estimated  frequency  and  time 
resolution  parameters  were  computed  from  the  F-T  plot  in  Figure  4.2  as  Af  «  36.62  KHz 
and  At  «  1.36xl0'5  s.  It  can  be  seen  that  the  signal  frequency  starts  and  remains  at 
/=  95  KHz  until  t  =  0.65  ms,  at  which  time  the  frequency  decreases  linearly  to  t  =  1  ms. 
At  t=\  ms  the  frequency  is  /=  45  KHz  which  remains  constant  until  the  end  of  the 
signal.  It  also  can  be  seen  that  the  peak  signal  power  is  higher  in  the  lower  frequency  of 
/=  45  KHz  and  for /=  95  KHz  the  average  power  is  approximately  -3.0  dB  compared  to 
the  maximum  signal  power.  By  comparison  with  analytic  signal  responses  shown  in 
Figure  4.1,  the  QMFB  frequency  and  time  averages  show  that  the  process,  with  some 
degradation  and  signal  processing  loses  due  to  the  computing  processing  and  the  instant 
changes  of  signal  frequency,  are  consistent. 
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Figure  4.2.  QMFB  Layer  #14  output  for  LFM  V&V  signal.  Average  PSD  plot  based  on 
Layer  #14  as  presented  and  average  time  plot  based  on  Layer  #10. 

4.2.2  Analytic  Szmajda  V&V  Signal 

The  normalized  time  and  frequency  responses  for  the  analytic  Szmajda  signal  are 
presented  in  Figure  4.3.  The  signal  was  resampled  and  zero-padded  prior  to  QMFB 
processing  according  to  the  values  shown  in  Table  4.2.  Related  to  the  time  response,  four 
different  signal  magnitudes  can  be  seen.  Besides  the  absence  of  a  signal  response  in  the 
t  =  25  s  to  t  =  37  s  interval,  no  other  noticeable  parameters  can  be  identified.  Related  to 
the  PSD  response,  it  can  be  seen  that  considering  the  W-sodB  bandwidth,  three  carrier 
frequencies  are  present,  including  fc  =  5  Hz,  f=  10  Hz  and  f:  =  40  Hz.  The  signal 
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component  located  at  fc  =  5  Hz  contains  the  higher  power  of  (P  =  0  dB)  compared  to  the 
equal  power  components  located  at  fc=  10  Hz  and  /,’  =  40  Hz  (P  =  -8  dB,  each). 


Frequency  (Hz) 


Figure  4.3  Szmajda  V&V  signal  time  domain  and  frequency  responses 


Table  4.2.  Szmajda  V&V  Signal  Parameters 


Input 

Bandwidth 

(Hz) 

Samp  Freq 
fs  (Hz) 

Duration 

(Sec) 

Number  of 
Samples 

Sample 

Rate 

Zero 

Padding 

Original 

4.00x1 04 

l.OOxlO2 

5-OOxlO1 

5000 

1 

N/A 

Resampled 

4.00x1 04 

3.00xl05 

1.12xl02 

15000000 

3000 

225 

After  pre-processing  the  Szmajda  signal  according  to  parameters  in  Table  4.2,  the 
signal  was  input  to  the  QMFB  process.  Representative  QMFB  results  for  Layer  #18  are 
presented  in  Figure  4.4  which  shows  the  discrete  frequency  responses  in  the  signal.  It  can 
be  seen  that  the  signal  started  at /=  5.0  KHz  and  increased  its  power  in  the  interval 
13.0  <  t  <  25.0  s.  For  the  interval  25.0  <t<  37.0  s  there  is  no  signal  present,  followed  by 
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the  interval  37.0  <t<  50.0  s  from  37  s  to  50  s  when  three  frequencies  are  present 
(/"=5.0KHz,  /=  10.0  KHz,  and  /=  40.0  KHz).  The  estimated  frequency  and  time 
resolution  parameters  were  computed  from  the  F-T  plot  in  Figure  4.4  as  Af  «  0.57  KHz 
and  At  &  0.88  s.  For  the  average  frequency  and  time  plots,  the  parameters  are 
Af  ~  0.57  KHz  and  At  «  13.65  ms,  respectively.  By  comparing  analytic  signal  responses 
in  Figure  4.3,  both  the  average  frequency  and  time  responses  from  the  QMFB  process, 
with  some  degradation  and  loses  due  to  the  computing  processing,  was  consistent. 


Figure  4.4.  QMFB  Layer  #18  output  for  Szmajda  V&V  signal.  Average  PSD  plot  based 
on  Layer  #18  as  presented  and  average  time  plot  based  on  Layer  #12. 
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After  processing  and  analyzing  QMFB  perfonnance  using  the  V&V  signals,  it 
was  concluded  that  the  QMFB  process  is  able  to  effectively  process  different  types  of 
signals  and  produce  outputs  matching  theoretical  expected  results.  However,  some  pre- 
QMFB  filtering  artifacts  appeared  due  to  the  effects  of  instantaneous  frequency  on  filter 
performance.  Beside  the  artifacts,  the  QMFB  results  presented  an  acceptable  response  in 
frequency  and  time  allocation  within  the  middle  layers. 

4.3  OFDM-Based  Signal  Performance 

The  two  OFDM-based  signals  described  in  Chapter  2  were  used  as  input  to  the 
QMFB  process  and  simulation  perfonned  according  to  the  methodology  explained  in 
Chapter  3.  Results  for  802.11a  Wi-Fi  and  802. 16e  WiMAX  OFDM-based  signals  are 
presented  and  QMFB  output  reliability  assessed  relative  to  input  signal  features.  Overlay 
plots  are  computed  to  compare  input  signal  and  QMFB  output  responses  for  single  burst 
and  integrated  burst  response  cases. 

4.3.1  Experimental  802.11a  Wi-Fi  Signal 

Experimental  802.11a  Wi-Fi  signal  assessment  was  performed  using  data 
collected  in  support  of  previous  work  detailed  in  [8].  Collections  were  made  using 
AFIT’s  RFSICS  with  the  Wi-Fi  devices  operating  in  an  anechoic  chamber  environment. 
According  to  the  structure  of  this  signal  described  in  Chapter  2,  the  802.11a  signal  is 
composed  of  two  distinct  regions.  The  preamble  region  is  used  for  network 
synchronization,  timing,  control,  etc.,  and  the  payload  region  is  used  for  transferring  user 
data.  Per  IEEE  standards  for  802.1  la  implementation  [9],  the  preamble  is  further  divided 
into  two  distinct  regions,  with  the  first  half  containing  10  short  OFDM  symbols  and  the 


35 


second  half  containing  2  long  OFDM  symbols.  The  10  short  OFDM  symbols  region  is 
selected  here  for  demonstration  with  a  single  burst  sent  to  the  QMFB  process  first  and 
then  an  integrated  collection  of  bursts  sent  to  the  QMFB  process. 

4. 3. 1.1  802.1  la  Wi-Fi  Preamble 

The  normalized  time  and  PSD  responses  for  802.1  la  Wi-Fi  preamble  are  shown 
in  Figure  4.5.  The  first  ten  symbols  of  the  preamble  were  isolated,  resampled  and  zero- 
padded  prior  to  QMFB  processing  according  to  the  values  shown  in  Table  4.3.  Related  to 
the  time  response,  the  signal  duration  is  approximate  T=  8  ps  and  it  can  be  seen  that  10 
peaks  are  very  noticeable.  Considering  the  W-wm  bandwidth,  the  PSD  response  shows 
twelve  distinct  frequency  components  that  match  the  802.11a  signal  structure  described 
in  Chapter  2. 


Figure  4.5  Time  and  PSD  responses  for  10  short  symbols  802.1  la  Wi-Fi  preamble. 
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Table  4.3.  802.1  la  10  short  symbols  set  as  input  to  QMFB 


Input 

Bandwidth 

(Hz) 

Samp  Freq 
fs  (Hz) 

Duration 

(Sec) 

Number  of 
Samples 

Sample 

Rate 

Zero 

Padding 

Original 

8.00xl06 

2.38xl07 

1.35xl0'4 

3200 

1 

N/A 

Resampled 

8.00xl06 

2.38xlOn 

1.41xl0'4 

32000000 

10000 

225 

The  QMFB  output  for  Layer  #19  is  shown  in  Figure  4.6.  The  allocations  of  twelve 
frequency  components  seen  in  Figure  4.5  are  present.  The  estimated  frequency  and  time 
resolution  parameters  were  computed  from  the  F-T  plot  in  Figure  4.6  as  Af  «  1 13.25 
KHz  and  At  ~  4.56  ps.  For  the  average  frequency  and  time  plots,  the  parameters  are 
Af  ~  56.6  KHz  and  At  ~  0.069  ps,  respectively.  Appendix  C.l  presents  a  table  of 
computed  Af  and  At  resolutions  as  function  of f  for  a  given  layer.  It  can  be  seen  that  the 
frequencies  are  allocated  within  the  expected  interval  (f=  1  MHz  to  /=  8  MHz).  The 
average  time  plot  was  computed  using  Layer  #14  due  to  the  F-T  trade-off  explained  in 
Chapter  2  with  better  time  resolution  obtained  in  lower  QMFB  layers.  The  expected  ten 
peaks  are  noticeable  in  the  average  time  plot  and  are  located  approximately  at  the  same 
expected  times.  However,  the  output  signal  magnitude  response  is  between  1.0  and 
2.0  dB  lower  when  compared  to  the  input  signal  in  Figure  4.5. 


37 


I 


■10 


••15 


•20 


-25 


'■30 


Figure  4.6  QMFB  Layer  #20  output  for  10  short  symbols  of  802. 11a  Wi-Fi  preamble. 
Average  PSD  based  on  Layer  #21  as  presented  and  average  time  based  on  Layer  #14. 


4.3. 1.2  802.11a  Wi-Fi  Single  Burst  Response 

The  normalized  time  and  PSD  responses  using  a  single  802.11a  Wi-Fi  burst  are 
shown  in  Figure  4.7.  The  signal  was  resampled  and  zero-padded  prior  to  QMFB 
processing  according  to  the  values  shown  in  Table  4.4.  Related  to  the  time  response,  the 
signal  duration  is  approximate  T=  11.6  ms  with  ten  visible  peaks  corresponding  to  the 
first  half  of  the  preamble.  It  can  be  seen  that  the  rest  of  the  signal  presents  a  noise  like 
behavior.  The  PSD  response  shows  that  the  signal’s  frequency  components,  considering 
the  W-3 odB  baseband  bandwidth,  the  dominant  frequency  responses  are  located  from /=  0 
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to /=  10  MHz.  These  initial  results  matched  the  signal  structure  described  in  Chapter  2. 


Table  4.4.  802. 11a  single  burst  signal  values  set  as  input  to  QMFB 


Input 

Bandwidth 

(Hz) 

Samp  Freq 
fs  (Hz) 

Duration 

(Sec) 

Number  of 
Samples 

Sample 

Rate 

Zero 

Padding 

Original 

9.00xl06 

2.38xl07 

1.24x1  O'4 

2945 

1 

N/A 

Resampled 

9.00xl06 

2.38xl010 

1.77xl0'4 

2945000 

1000 

222 

The  QMFB  output  for  Layer  #16  is  shown  in  Figure  4.8.  The  allocations  of 
multiple  frequency  components  seen  in  Figure  4.7  are  present.  The  estimated  frequency 
and  time  resolution  parameters  were  computed  from  the  F-T  plot  in  Figure  4.8  as 
Af  *  181.12  KHz  and  At « 2.8  ps.  For  the  average  frequency  and  time  plots,  the 
parameters  are  Af  «  1 1.3  KHz  and  At  *21.56  ns,  respectively.  Appendix  C.2  presents  a 
table  of  computed  Af  and  At  resolutions  as  function  of  f  for  a  given  layer.  It  can  be  seen 
that  frequencies  are  allocated  within  the  expected  interval  {f—  1  MHz  to  /=  10  MHz). 
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The  first  half  of  the  preamble  response  is  still  noticeable  and  is  located  approximately  at 
the  same  expected  times  compared  to  the  input  shown  in  Figure  4.7. 


According  to  the  802.11a  signal  parameters  described  in  Chapter  2,  there  are 
some  signal  features  that  can  be  extracted  through  qualitative  visual  assessment  of  the 
QMFB  output  shown  in  Figure  4.8: 

•  The  10  short  symbols  shown  in  the  left-hand  red  dashed  rectangle 

•  The  guard  Interval  shown  in  the  right-hand  red  dashed  rectangle 


Figure  4.8  QMFB  Layer  #16  output  for  single  802. 1  la  Wi-Fi  burst.  Average  PSD  based 
on  Layer  #20  as  presented  and  average  time  based  on  Layer  #9. 
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4.3.1 .3  Wi-Fi  Integrated  Burst  Response 

The  normalized  time  and  PSD  responses  for  integrated  802.11a  bursts  are 
shown  in  Figure  4.9.  In  this  case,  a  total  of  As  =  500  802.11a  burst  responses  were 
integrated  to  create  a  new  input  to  the  QMFB  process.  The  signal  was  resampled  and 
zero-padded  prior  to  QMFB  processing  according  to  the  values  shown  in  Table  4.5. 
Related  to  the  time  response,  the  signal  duration  is  approximate  T=  11.6  ms  with  ten 
visible  peaks  corresponding  to  the  first  half  of  the  preamble  and  the  rest  of  the  signal 
presents  a  noise  like  behavior.  The  PSD  response  shows  that  the  signal’s  frequency 
components,  considering  the  IV.  mm  baseband  bandwidth,  the  dominant  frequency 
responses  are  located  from /=  0  to  /=  10  MHz.  These  results  match  the  signal  structure 
described  in  Chapter  2.  It  also  can  be  seen  that  the  first  ten  symbols  show  a  uniform 
magnitude,  with  the  peak  located  at  approximately  at  t  =  2.4  ps  corresponding  to  the 
guard  interval  described  in  Chapter  2.  The  part  of  the  signal  corresponding  to  the  payload 
(user  data)  shows  a  magnitude  reduction  due  to  randomness  of  the  symbol  assignment. 
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Figure  4.9  Time  and  PSD  responses  for  Nb  =  500  integrated  bursts  802. 11a  Wi-Fi. 


Table  4.5.  802. 1  la  signal  values  for  Nb  =  500  integrated  bursts  as  input  to  QMFB 


Input 

Bandwidth 

(Hz) 

Samp  Freq 
fs  (Hz) 

Duration 

(Sec) 

Number  of 
Samples 

Sample 

Rate 

Zero 

Padding 

Original 

9.00xl06 

2.38xl07 

1.24x1  O'4 

2945 

1 

N/A 

Resampled 

9.00xl06 

2.38xl010 

1.77x1 0'4 

2945000 

1000 

222 

The  QMFB  output  for  Layer  #16  is  shown  in  Figure  4.10for  integration  of 
Nb  =  500  bursts.  The  allocation  of  multiple  frequencies  components  seen  in  Figure  4.9  is 
present.  The  estimated  frequency  and  time  resolution  parameters  were  computed  from  the 
F-T  plot  in  Figure  4.10  as  Af  «  181.12  KHz  and  At «  2.8  ps.  For  the  average  frequency 
and  time  plots,  the  parameters  are  Af  «  22.6  KHz  and  At  «  0. 1 73  ns,  respectively. 
Appendix  C.2  presents  a  table  of  computed  Af  and  At  resolutions  as  function  of  f  for  a 
given  layer.  The  first  half  of  the  preamble  response  is  still  noticeable  and  is  located 
approximately  at  the  same  expected  times  compared  to  the  input  shown  in  Figure  4.9. 
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According  to  the  802.11a  signal  parameters  described  in  Chapter  2,  there  are 
some  signal  features  that  can  be  extracted  through  qualitative  visual  assessment  of  the 


QMFB  output  shown  in  Figure  4.10: 

•  The  10  short  symbols  shown  in  the  left-hand  red  dashed  rectangle 

•  Lower  correlation  in  the  payload  (user  data)  region  in  the  right-hand  red 
dashed  rectangle 

•  Average  signal  power  below  -5.0  dB  for  the  payload  (user  data)  region 
highlighted  by  a  red  arrow  in  the  average  time  plot. 


Figure  4.10  QMFB  Layer  #16  output  for  Nb  =  500  integrated  Wi-Fi  bursts.  Average 
PSD  plot  based  on  Layer  #19  as  presented  and  average  time  plot  based  on  Layer  #12. 
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4.3. 1.4  802.1  la  Wi-Fi  Preamble:  Single  vs.  Integrated  Response 

To  assess  the  QMFB  response  relative  to  signal  input  features,  average  time  and 
frequency  plots  are  provided  in  Figure  4.1 1  to  Figure  4.14  for  single  burst  and  integrated 
burst  QMFP  processing.  Related  to  the  overlaid  time  responses  in  Figure  4.11  and  Figure 
4.12,  it  can  be  seen  that  the  QMFB  output  envelope  matches  the  input  time  response  for 
both  single  and  integrated  burst  cases.  Therefore,  the  F-T  plots  computed  during  the 
process  are  reliable  for  revealing  802.11a  Wi-Fi  signal  characteristics.  In  the  payload 
(user  data)  region  of  the  average  time  response  ( t  >  200  ps),  the  power  reduction  for  burst 
integration  is  evident  given  the  random  signal  structure  in  this  region.  Related  to  the 
average  PSD  response  presented  in  Figure  4.13  and  Figure  4.14,  the  burst  integration 
resulted  in  gain  of  approximate  G  «  3.0  dB  when  compared  to  the  input  signal  or  single 
burst  QMFB  responses. 
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Figure  4.12  Average  time  responses  for  802. 11a  Wi-Fi  signal  expanded  region  for 

0  <  t<  200  ps. 


Figure  4.13  Average  PSD  responses  for  802.1  la  Wi-Fi  signal. 
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Figure  4. 14  Average  PSD  responses  for  802. 11a  Wi-Fi  signal  expanded  region  for 

0  </<  3  MHz. 


4.3.2  Experimental  802. 16e  WiMAX  Signal 

Experimental  802. 16e  WiMAX  signal  assessment  was  performed  using  data 
collected  in  support  of  previous  work  detailed  in  [0,  11].  The  collections  were  obtained 
using  AFIT’s  RFSICS  with  Alvarion  BreezeMAX  5000  Mobile  Subscriber  (MS)  devices 
operating  in  a  typical  office  environment  [13].  According  to  the  structure  of  this  signal 
described  in  Chapter  2,  and  experimental  observations  noted  in  [0,  11],  the  analysis  was 
divided  in  two  parts.  The  first  part  only  considers  WiMAX  range-only  burst  responses 
and  the  second  part  considers  only  WiMAX  data-only  burst  responses.  For  each  of  these 
cases,  QMFB  processing  is  conducted  using  single  burst  and  integrated  burst  responses. 
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4.3.2. 1  WiMAX  Range-Only  Burst:  Single  Response 

The  normalized  time  and  PSD  responses  using  a  single  range-only  WiMAX 
burst  are  shown  in  Figure  4.15.  The  signal  was  resampled  and  zero-padded  prior  to 
QMFB  processing  according  to  the  values  shown  in  Table  4.6.  Related  to  the  time 
response,  the  signal  duration  is  approximately  7=0.4  ms  and,  considering  the  W~. sodB 
bandwidth,  the  PSD  response  clearly  shows  some  frequency  components  located  in  three 
spectral  regions,  including:  0  </<  0.4  MHz,  0.6  </ <1.0  MHz  and  1.6  </ <  2.0  MHz. 


Frequency  [Hz] 


x  106 


Figure  4.15.  Time  and  PSD  responses  for  single  802. 16e  WiMAX  range-only  burst. 


Table  4.6.  802. 16e  WiMAX  range-only  single  burst  parameters 


Input 

Bandwidth 

(Hz) 

Samp  Freq 
fs  (Hz) 

Duration 

(Sec) 

Number  of 
Samples 

Sample 

Rate 

Zero 

Padding 

Original 

2.00xl06 

1.19xl07 

3.81xl0'4 

4521 

1 

N/A 

Resampled 

2.00xl06 

3.56xl010 

9.42x1  O'4 

13563000 

3000 

225 
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The  QMFB  output  for  Layer  #18  is  shown  in  Figure  4.16  for  single  802. 16e 
WiMAX  range-only  burst.  The  allocation  of  six  frequencies  components  seen  in  Figure 
4.15  is  present.  The  estimated  frequency  and  time  resolution  parameters  were  computed 
from  the  F-T  plot  in  Figure  4.16  as  Af  «  67.9  KHz  and  At  &  7.4  ps.  For  the  average 
frequency  and  time  plots,  the  parameters  are  Af  «  16.9  KHz  and  At  ~  0. 1  1 5  ps, 
respectively.  Appendix  C.3  presents  a  table  of  computed  Af  and  At  resolutions  as  function 
of  fs  for  a  given  layer.  It  can  be  seen  that  frequencies  are  allocated  within  the  three 
expected  intervals  (0</<0.4,  0.6  </<  1.0,  and  1.6</<2.0  MHz),  are  distributed  in 
three  pairs,  and  span  the  bandwidth  shown  in  Figure  4.15.  Within  each  pair  of 
frequencies  a  transition  is  seen  between  a  certain  numbers  of  transmitted  symbols. 
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Figure  4.16  QMFB  Layer  #18  output  for  single  802. 16e  WiMAX  range-only  burst. 
Average  PSD  plot  based  on  Layer  #20  and  average  time  plot  based  on  Layer  #12. 


4. 3. 2.2  WiMAX  Range-Only  Burst:  Integrated  Response 

The  normalized  time  and  PSD  responses  using  integrated  range-only  WiMAX 
bursts  are  shown  in  Figure  4.17.  In  this  case,  a  total  of  Nb  —  1400  802. 16e  burst 
responses  were  integrated  to  create  a  new  input  to  the  QMFB  process.  The  range-only 
integrated  burst  signal  was  resampled  and  zero-padded  prior  to  QMFB  processing 
according  to  the  same  values  previously  shown  in  Table  4.6.  Related  to  the  time 
response,  the  signal  duration  is  approximately  r=0.4  ms  and,  considering  the  W- aodB 
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bandwidth,  the  PSD  response  contains  the  same  frequency  components  as  the  single  burst 
response  (0  <  /<0.4  MHz,  0.6</<1.0  MHz  and  1.6</<2.0  MHz)  plus  some 
additional  frequency  components  below  the  P  =  -20  dB  level. 


Figure  4. 17.  Time  and  PSD  responses  for  NB  =  1400  integrated  802. 16e  WiMAX  range- 

only  bursts. 


The  QMFB  output  for  Layer  #18  is  shown  in  Figure  4.18  for  Nb  =1400 
integrated  802. 16e  WiMAX  range-only  bursts.  The  allocation  of  six  frequencies 
components  seen  in  Figure  4.17  is  present.  The  estimated  frequency  and  time  resolution 


parameters  were  computed  from  the  F-T  plot  in  Figure  4.18  as  Af  «  67.9  KHz  and 


At  ~  7.4  us.  For  the  average  frequency  and  time  plots,  the  parameters  are  Af  «  16.9  KHz 


and  At»  0.115  lis,  respectively.  Appendix  C.3  presents  a  table  of  computed  Af  and  At 
resolutions  as  function  of f  for  a  given  layer.  It  can  be  seen  that  frequencies  are  allocated 
within  the  three  expected  intervals  (0  </<  0.4,  0.6  </<  1.0,  and  1.6  </<  2.0  MHz),  are 
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distributed  in  three  pairs,  and  span  the  bandwidth  shown  in  Figure  4.17.  In  this  case,  the 
middle  pair  of  frequencies  presented  the  higher  power  correlation  compared  to  the  other 
frequency  components. 


Figure  4.18  QMFB  Layer  #18  output  for  Nb  =  1400  integrated  802. 16e  WiMAX  range- 
only  bursts.  Average  PSD  based  on  Layer  #20  and  average  time  based  on  Layer  #12. 

4. 3.2. 3  WiMAX  Range-Only  Burst:  Single  vs.  Integrated  Response 

To  assess  the  QMFB  response  relative  to  input  signal  features,  average  time  and 
frequency  plots  are  provided  in  Figure  4.19  to  Figure  4.22  for  single  burst  and  integrated 
burst  QMFP  processing.  Related  to  the  overlay  time  response  in  Figure  4.19  and  Figure 
4.20,  the  integrated  burst  QMFB  output  show  less  energy  in  the  time  response  due  to 
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integration  of  random  symbol  assignment  within  and  across  every  integrated  burst. 
Related  to  overlay  PSD  responses  in  Figure  4.21  and  Figure  4.22,  it  can  be  seen  that  the 
JF-3(mb  bandwidth  remains  constant  and  that  burst  integration  resulted  in  “thinner” 
responses  for  the  six  frequency  components. 


■Time  [s]  x  1Q-4 

Figure  4.19  Average  time  responses  for  802. 16e  WiMAX  range-only  burst. 
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Figure  4.20  Average  time  responses  for  802. 16e  WiMAX  range-only  burst.  Expanded 

region  for  0  <  t  <  200  pis. 


1 

1 

. 

Input 

Single  Burst 
Burst  Integration 

0.2  0.4  0.6  0.8  1  1.2  1.4  1.6  1.8  2 

Frequency  [Hz]  x  106 

Figure  4.2 1  Average  PSD  responses  for  802. 16e  WiMAX  range-only  burst. 
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Figure  4.22  Average  PSD  responses  for  802. 16e  WiMAX  range-only  burst.  Expanded 

region  for  0  </<  0.5  MHz. 

4. 3. 2.4  802. 16e  WiMAX  Data-Onlv  Burst:  Single  Response 

The  normalized  time  and  PSD  responses  using  a  single  data-only  802. 16e 
WiMAX  burst  is  shown  in  Figure  4.23.  The  signal  was  resampled  and  zero-padded  prior 
to  QMFB  processing  according  to  the  values  shown  in  Table  4.7.  Related  to  the  time 
response,  the  signal  duration  is  approximate  7=  1.6  ms  and  the  noise  like  behavior  of  the 
signal  is  seen  in  both  time  and  frequency  responses. 
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Figure  4.23  Time  and  PSD  responses  for  single  802. 16e  WiMAX  data-only  burst. 


Table  4.7.  802.  16e  WiMAX  data-only  single  burst  parameters 


Input 

Bandwidth 

(Hz) 

Samp  Freq 
fs  (Hz) 

Duration 

(Sec) 

Number  of 
Samples 

Sample 

Rate 

Zero 

Padding 

Original 

2.50xl06 

1.19xl07 

1.61xl0'3 

19125 

1 

N/A 

Resampled 

2.50xl06 

1.19xl010 

2. 83x1 0'3 

19125000 

1000 

225 

The  QMFB  output  for  Layer  #16  is  shown  in  Figure  4.24  for  a  single  802. 16e 
WiMAX  range-only  burst.  The  multiple  frequency  allocations  shown  in  Figure  4.23  are 
present.  The  estimated  frequency  and  time  resolution  parameters  were  computed  from 
the  F-T  plot  in  Figure  4.24  as  Af  «  90.6  KHz  and  At  &  5.53  ps.  For  the  average 
frequency  and  time  plots,  the  parameters  are  Af  «  5.66  KHz  and  At  ~  0.345  us, 
respectively.  Appendix  C.3  presents  a  table  of  computed  Af  and  At  resolutions  as  function 
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of fs  for  a  given  layer.  Other  than  bandwidth  (W. 30dB  ~  2.25  MHz)  and  signal  duration 


(T  x  1.6  ms),  no  signal  structure  can  be  extracted  from  the  F-T  plot. 


x  lO'3 


Figure  4.24  QMFB  Layer  #16  output  for  single  802. 16e  WiMAX  data-only  burst. 
Average  PSD  based  on  Layer  #20  and  average  time  based  on  Layer  #12. 

4. 3. 2. 5  802. 16e  WiMAX  Data-Only  Burst:  Integrated  Response 

The  normalized  time  and  PSD  responses  for  integrated  data-only  WiMAX 

bursts  are  shown  in  Figure  4.25.  The  signal  was  resampled  and  zero-padded  prior  to 

QMFB  processing  according  to  the  values  previously  shown  in  Table  4.7.  Related  to  the 

time  response,  the  signal  duration  is  approximate  1.6  ms  and  it  can  be  seen  the  noise  like 

behavior  of  the  signal  in  both  time  and  frequency  responses. 
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Frequency  [Hz]  x  1Q« 

Figure  4.25  Time  and  PSD  responses  for  Nb  =  640  integrated  802. 16e  WiMAX  data- 

only  bursts. 

The  QMFB  output  for  Layer  #16  is  shown  in  Figure  4.26  for  Nb  =  640  integrated 
802. 16e  WiMAX  data-only  bursts.  The  allocation  of  multiple  frequencies  components 
seen  in  Figure  4.25  is  present.  The  estimated  frequency  and  time  resolution  parameters 
were  computed  from  the  F-T  plot  in  Figure  4.26  as  Af  «  90.6  KHz  and  At  ~  5.53  ps.  For 
the  average  frequency  and  time  plots,  the  parameters  are  Af  «  5.66  KHz  and  At  ~  0.345 
ps,  respectively.  Appendix  C.3  presents  a  table  of  computed  Af  and  At  resolutions  as 
function  of  f  for  a  given  layer.  Other  than  bandwidth  (W. 30dB  ~  2.25  MHz)  and  signal 
duration  (T«  1.6  ms),  no  signal  structure  can  be  extracted  from  the  F-T  plot. 
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Figure  4.26  QMFB  Layer  #16  output  for  Nb  =  640  integrated  802. 16e  WiMAX  data- 
only  bursts.  Average  PSD  based  on  Layer  #20  and  average  time  based  on  Layer  #12. 

4. 3. 2. 6  802. 16e  WiMAX  Range-Only  Burst:  Single  vs.  Integrated  Response 

To  assess  the  QMFB  response  relative  to  input  signal  features,  average  time  and 
frequency  plots  are  provided  in  Figure  4.27  to  Figure  4.30  for  single  burst  and  integrated 
burst  QMFP  processing.  Related  to  the  overlaid  time  responses  in  Figure  4.27  and  Figure 
4.28  the  integrated  burst  QMFB  output  shows  less  energy  in  the  time  response  due  to 
integration  of  random  symbol  assignment  within  and  across  every  integrated  burst. 
Related  to  the  PSD  responses  in  Figure  4.29  and  Figure  4.30,  it  can  be  seen  that  the 
bandwidth  is  constant  along  the  burst  integration. 
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Figure  4.27  Average  time  responses  for  802. 16e  WiMAX  data-only  burst 


o 


Figure  4.28  Average  time  responses  for  802. 16e  WiMAX  data-only  burst.  Expanded 

region  for  0  <  t  <  200  jus  (Bottom). 
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Figure  4.29  Average  PSD  responses  for  802. 16e  WiMAX  data-only  burst. 


Figure  4.30  Average  PSD  responses  for  802. 16e  WiMAX  data-only  burst.  Expanded 

region  for  0  </<  1 .0  MFIz. 


60 


4.4  Summary 

The  results  presented  in  this  chapter  enabled  qualitative  visual  assessment  of 
QMFB  performance  using  each  signal  of  interest.  Related  to  the  LFM  and  Szmajda  V&V 
signals,  it  was  seen  that  the  QMFB  output  offers  reliable  estimation  of  frequency  and 
time  characteristics.  Once  the  baseline  was  defined,  verified  and  validated,  experimental 
OFDM-based  signals  were  introduced  and  QMFB  perfonnance  assessed.  MATLAB® 
simulation  results  and  qualitative  visual  analysis  were  sufficient  for  identifying  some 
signal  features  such  as  signal  duration  (Ts)  and  bandwidth  ( W.  vidis)  while  highlighting 
relative  signal  structure. 

Related  to  the  analyzed  802.1  la  Wi-Fi  signal,  the  preamble  structure  revealed  by 
QMFB  processing  is  consistent  with  the  theoretical  description  in  Chapter  2.  Using  burst 
integration,  the  effect  of  random  symbol  assignment  during  the  payload  (user  data)  region 
became  evident  as  lower  magnitude  responses  in  Figure  4.10.  The  average  time  and 
average  frequency  plot  showed  that  the  QMFB  output  envelope  followed  the  original 
signal  behavior,  but  average  loses  of  approximate  1.0  dB  for  burst  integration  and  2.0  dB 
for  single  burst  were  observed. 

Related  to  the  analyzed  802. 16e  WiMAX  signal,  the  results  are  consistent  with 
the  framework  defined  in  Chapter  2  and  the  QMFB  output  for  both  the  ranging-on  and 
data-only  cases  enabled  qualitative  visual  assessment  and  provided  some  exploitation 
benefit.  Some  key  signal  features  such  as  bandwidth,  frequency  and  time  allocation  were 
extractable,  and  after  further  layer  analysis  approximations  to  symbol  duration  and 
channel  bandwidth  were  available. 
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CHAPTER  5. 


Summary,  Conclusions  and  Recommendations 

This  chapter  presents  a  summary  and  conclusions  of  research  findings  along  with 
some  recommendations  for  future  work  that  have  been  extracted  as  result  of  this  research 
effort.  The  chapter  is  divided  in  three  subsections:  1)  Section  5.1  presents  a  summary 
review  related  to  the  topics  of  interest  covered  during  the  research  with  some  key 
exploitation  concepts  highlighted  relative  to  the  overall  process;  2)  Section  5.2  presents 
conclusions  extracted  from  results  after  applying  the  adopted  methodology  in  Chapter  3 
and  the  baseline  definitions  given  in  chapter  2;  and,  3)  Section  5.3  is  aimed  at  identifying 
other  applications  for  QMFB  exploitation  providing  recommendations  for  future  work. 

5.1  Summary 

Chapter  1  presented  the  operational  and  technical  motivation  along  with  the 
research  objectives.  Related  to  the  operational  motivation,  the  approach  was  to  provide  a 
qualitative  visual  assessment  process  for  a  given  signal  from  an  operator’s  perspective. 
By  perfonning  QMFB  processing  with  a  given  signal  it  is  possible  to  obtain  frequency¬ 
time  (F-T)  plots  [1,  3,  4]  having  features  that  represent  signal  features  such  as  bandwidth 
( W ),  center  frequency  (f  ),  signal  duration  (Ts),  modulation  type  (AM,  FM,  BPSK,  QAM, 
etc),  frequency  content  and  time  allocation.  Because  of  its  unique  structure,  two  widely 
used  signals  based  on  Orthogonal  Frequency  Division  Multiplexing  (OFDM)  were 
chosen  as  signals  of  interest  for  demonstration.  The  operational  motivation  relies  on 
optimizing  the  QMFB  process  for  the  two  signals  of  interest  having  defined  structure 
according  to  [9,  12]. 
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Chapter  2  presented  the  necessary  technical  background  aimed  at  describing  the 
topics  of  interest  related  to  the  operational  and  technical  motivation.  The  general 
implementation  of  the  QMFB  process  was  described  and  implementation  parameters  and 
assumptions  were  made.  The  basic  structure  of  OFDM  signals  was  introduced  next 
related  to  the  physical  layer  perspective  of  802.011a  Wi-Fi  and  802. 16e  WiMAX  frame 
structures. 

Chapter  3  presented  the  adopted  methodology  aimed  at  satisfying  the  defined 
research  objectives  using  the  technical  details  described  in  Chapter  2.  An  overall  process 
overview  was  presented  and  the  Verification  and  Validation  (V&V)  signals  were 
described  [2,  7].  The  802.011a  Wi-Fi  and  802. 16e  WiMAX  experimental  signals 
parameters  were  addressed  [8,  0,  11,  13]  and  the  initial  QMFB  configuration  was 
established.  The  effect  of  resampling  and  zero-padding  was  described  and  measurable 
outputs  were  defined.  Finally,  the  graphical  presentation  fonnat  for  QMFB  results  was 
introduced. 

Chapter  4  presented  computed  simulation  results  that  were  obtained  after  applying 
the  defined  methodology  to  each  signal  of  interest.  Initial  time  domain  and  frequency 
domain  responses  were  presented  for  each  input  signal  along  with  the  initial  and 
computed  resampled  parameters  for  each  case.  Results  for  selected  QMFB  outputs  were 
presented  according  to  the  process  introduced  in  Chapter  3  using  2D  T-F  QMFB  plots 
and  ID  average  frequency  and  average  time  plots.  Overlay  plots  of  the  ID  responses 
were  computed  and  the  input  signal  and  output  QMFB  responses  compared  using  a  single 
burst  and  multiple  integrated  bursts.  Time  window  resolution  (At)  and  frequency  window 
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resolution  (Af)  parameters  were  extracted  from  each  computed  QMFB  layer  output. 

5.2  Conclusions 

Computed  QMFB  outputs  were  generated  for  all  interest  signals  according  to  the 
adopted  methodology.  The  QMFB  process  was  implemented  and  decomposed 
successfully  according  to  the  defined  measurable  outputs.  The  following  conclusions  are 
presented  based  on  the  defined  measurable  outputs  presented  in  Chapter  3: 

1 .  Layer  Generation  Minimization 

Because  higher  number  QMFB  layers  did  not  improve  qualitative  analysis, 
or  did  not  give  any  additional  infonnation  due  to  coarser  At  resolution 
when  compared  with  lower  numbered  layers,  the  last  three  QMFB  layers 
were  not  computed.  Related  to  complete  layer  processing,  this  reduced 
overall  computation  time  by  approximate  78%  when  Nq  =  22  total  layers 
were  available  (Ns  =  2 ~  total  input  signal  samples)  and  83%  when  Nq  =  25 
total  layers  were  available  ( Ns  =  2 "  total  input  signal  samples). 
Representative  average  computing  times  are  shown  in  Table  5.1 


Table  5.1  Average  Layer  Computing  Time 


Layer 

# 

C.T. 

(min) 

Layer 

# 

C.T. 

(min) 

Layer 

# 

C.T. 

(min) 

Layer 

# 

C.T. 

(min) 

Layer 

# 

C.T. 

(min) 

1 

6 

6 

21 

11 

36 

16 

54 

21 

138 

2 

9 

7 

24 

12 

39 

17 

59 

22 

213 

3 

12 

8 

27 

13 

42 

18 

66 

23 

360 

4 

15 

9 

30 

14 

46 

19 

78 

24 

655 

5 

18 

10 

33 

15 

49 

20 

99 

25 

1242 
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2 .  Output  Matrix  Qn  Formatting 

After  computing  this  measurable  output  data  corresponding  to  At  and  Af 
sample  frequency,  signal  amplitude  (real  or  complex)  or  magnitude  values 
were  successfully  save  for  further  analysis  and  the  results  were  consistent 
with  studied  input  signals.  As  result  of  a  given  Qn  matrix  a  2D  F-T  plot 
was  computed  which  permitted  to  visually  assess  input  signal 
characteristics  such  as  bandwidth  ( W  -sodisf  signal  duration  (Ts),  time  and 
frequency  allocation. 

3.  Time  Resolution  (At)  Analysis 

This  output  was  used  to  ensure  consistency  between  the  average  QMFB 
time  vector  extracted  from  a  given  QN  matrix  and  the  input  signal  time 
response.  This  was  done  using  a  ID  overlay  plot  of  the  averaged  QMFB 
output  response  and  the  input  signal’s  time  response. 

4.  Frequency  Resolution  (Af)  Analysis 

This  output  was  used  to  ensure  consistency  between  the  average  QMFB 
frequency  vector  extracted  from  a  given  QN  matrix  and  the  input  signal 
frequency  domain  (PSD)  response.  This  was  done  using  a  ID  overlay  plot 
of  the  averaged  QMFB  output  response  and  the  input  signal’s  PSD 
response. 

After  evaluating  the  overall  performance  of  the  QMFB  process  with  the  signals  of 
interest,  it  can  be  concluded  that  it  is  possible  to  perfonn  reliable  qualitative  visual 
assessment  of  OFDM-based  signals.  Using  this  passive  method  some  distinctive  signal 
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parameters  can  be  extracted  such  as  bandwidth  (W),  center  frequency  (/■),  signal  duration 
(Ts),  frequency  content  and  time  allocation. 

Related  to  the  802.11a  Wi-Fi  signal,  the  10  short  preamble  symbols  were 
successfully  computed  and  results  of  the  QMFB  output  matched  the  standard  definition 
given  in  [9],  The  guard  interval  between  the  preamble  and  the  payload  was  also 
indentified  using  both  single  and  integrated  burst  inputs.  After  perfonning  burst 
integration,  the  QMFB  output  visually  matched  the  expected  results  using  2D  F-T  plots, 
with  lower  energy  observed  in  the  data  payload  region  due  to  random  symbol  assignment. 

Related  to  the  802. 16e  WiMAX  signal,  the  single  burst  range-only  and  data-only 
responses  were  successfully  processed  and  matched  the  standard  definition  given  in  [12]. 
Even  though  presenting  noise-like  signal  behavior  in  the  2D  F-T  plots,  QMFB  responses 
for  both  cased  were  sufficient  to  note  differences  in  frequency  and  time  allocations. 
Responses  using  integrated  burst  inputs  were  similar  as  those  computed  for  the  802.11a 
Wi-Fi  signal.  Key  signal  features  of  center  frequency  (fc),  frequency  content  and  time 
allocation  were  easily  extractable.  With  further  layer  analysis,  approximations  to  symbol 
duration  (Yv)  and  bandwidth  ( IV)  were  available  as  well. 

The  main  research  objective  was  met  and  the  QMFB  process  presents  an 
acceptable  response  for  extracting  some  key  features  from  802.011a  Wi-Fi  and  802. 16e 
WiMAX  signals.  This  empirical  approach  pennitted  exploitation  of  OFDM-based  signals 
using  qualitative  visual  assessment  such  as  may  benefit  human  analysts. 

5.3  Recommendations  for  Future  Research 

Further  exploitation  of  the  passive  QMFB  process  should  be  considered  in  future 
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work  based  on  different  QMFB  modification  and  decompositions.  Benefit  may  be 
realized  using  QMFB  processing  aimed  to  accomplish  automatic  signal  detection, 
channel  estimation  or  signal  fingerprinting  analysis.  Some  specific  future  research 
activities  could  include: 

1 .  OFDM-Based  Signal  Detection  Using  Correlation  Methods 

Create  a  signal  data  base  according  to  IEEE  standard  definitions  for  specific 
signals.  Simulated  data  base  features  could  be  correlated  with  experimental 
QMFB  exploited  features  with  a  goal  of  finding  distinctive  characteristics. 

2.  RF  Signal  Fingerprinting  Using  QMFB  Features 

Assess  the  use  of  experimental  2D  QMFB  (complex  or  real)  outputs  to 
generate  features  for  discriminating  between  signals  from  different  devices 
using  RF  fingerprints.  This  process  has  been  previously  demonstrated  using 
2D  Gabor  Transfonn  (GT)  features  [0,  11]  and  could  be  considered  for 
expanding  the  applicability  of  QMFB  processing. 

3.  Compare  2D  QMFB  Outputs  with  Other  2D  Transform  Outputs 
Improved  signal  exploitation  may  be  possible  using  different  2D  transforms, 
e.g.,  Gabor-based  [0,  11].  Performance  and  computation  time  could  be 
compared  with  QMFB  processing  as  implemented  here. 

4.  Consider  G-M  Transformed  QMFB  Processing 

Detection  and  estimation  may  be  improved  using  QMFB  outputs  and 
transformation  methods  investigated  previously  by  [4]  which  used  the 
Grohholz-Mims  (G-M)  2D  transfonnation  methods  for  non-cooperative 
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multi-channel  detection  [5,  6].  These  methods  have  not  been  applied  to 
OFDM-based  signals  such  as  done  here  with  QMFB  processing  and  results 
in  [4]  suggests  that  G-M  transfonned  QMFB  processing  may  be  well-suited 
for  additional  investigation. 
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APPENDIX  A  -  Simulation  Code 


This  appendix  presents  the  implemented  MATLAB®  code  according  to  the 
adopted  methodology  and  the  defined  measurable  outputs. 

A.  1  Time  and  Frequency  Input  Signal  Responses 

%  Original  file:  SigPSD_V2  Function 
%  Created  by:  Dr.  Michael  A.  Temple 

%  Modified  by:  Felipe  E.  Garrido,  CAP.  Chilean  Air  Force 

%  Last  update:  02/16/2012 


%  ======================  SigPSD_V3  Function 


function  [ SigPsd, Fscale, SigPow]  = 

SigPSD  V3 ( Sigln, Fsamp, SAve, NAve, SNorm, SPlot , t  lim, f  lim, dB  lim) 


Calculates  Row-Wise  Power  Spectral  Density  (PSD)  as  the  Magnitude 
of  Fourier  Transform  Coefficients  Squared.  Resultant  PSD  (SigPsd) 
and  Average  Power  Estimate  (SigPow)  are  output. 


Created:  Apr  06,  MAT 

Modified:  19  Apr  11,  MAT 
Modified:  7  Sep  11,  MAT  . 
Modified:  14  Dec  11,  MAT 


Enable  Matrix  Input/Output 
Minor  Change  to  Plotting  Output 
Save  time  and  frequency  responses  as 
vector . 


NOTE:  Outputs  Verif ied/Consistent  with  Matlab's  'PSD'  Function  for 

BOTH  Complex  and  Real-Valued  Input  Signal  (10  Dec  08) . 


Inputs : 
Sigln 
Fsamp 
SAve 

NAve 

SNorm 

SPlot 


t  lim 
f_lim 
dB  lim 


Nrow  x  NCol  Matrix:  NCol  =  Time  Domain  Signal  Samples 
Sample  Frequency  (Samples/Sec) 

PSD  Averaging  Control  Variable  (Sliding  Window  Average) 

(1  =  Average  Applied,  0  =  No  Average  Applied) 

Number  of  Points  Averaged  Across 

(Odd  Number  Required  if  Save  =  1,  Else,  Doesn't  Matter) 
Normalization  Control  Variable 
(1  =  Normalized,  0  =  UnNormalized) 

Plot  Control  Variable  (1  =  Plot,  0  =  No  Plot) .  Two  plots 
are  generated  for  SPlot=l :  A)  Time  Domain  Waveform 
(amplitude  for  Sigln  Real  and  Magnitude  for  Sigln  Complex) , 
B)  Sigln  PSD  (spanning  -Fsamp/2  <  f  <  Fsamp/2  for  Sigln 
Real  and  0  <  f  <  Fsamp  for  Sigln  Complex) . 

Time  window  =  to  plot  time  response  [0,t_lim] 

Frequency  window  =  to  plot  frequency  response  [0,f_lim] 
Magnitude  window  =  Defines  min  dB  level  to  plot  [dB  lim, 0] 


Outputs : 

SigPsd  PSD  of  Input  Signal  Calc  as  | Fourier{ Sigln} | A2 

Fscale  Resultant  Frequency  Scale  for  Output  PSD 

SigPow  Estimated  Power  Estimate  Calc  as  Sum  of  PSD  Components 
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function  [SigPsd, Fscale, SigPow]  = 

SigPSD_V3 ( Si gin, Fsamp, SAve, NAve, SNorm, SPlot , t_lim, f_lim, dB_lim) 

[Nrowln, NcolIn]=size (Sigln) ; 

if  NcolIn==l  %  Column  Vector  Input  Transpose  &  Process  as  Row  Vector 

SigIn=SigIn. ' ; 

[Nrow, Ncol]=size (Sigln) ; 

else 

Nrow=NrowIn; 

Ncol=NcolIn; 

end 

SigVal  =  isreal (Sigln (1, 1) ) ;  %  SigVal=l  (Real)  or  SigVal=0  (Complex) 

SigFFT= (fft (real (Sigln. ' ) ) . ' /Ncol) ; 

SigPsd=abs  (SigFFT)  .''2;  %  Non-Shifted  PSD 

SigPow=sum (SigPsd. ' ) ; 

if(SAve==l)  %  Apply  Row-Wise  Averaging/Smoothing  to  PSD 

NTest  =  mod(NAve,2);  %  Check/Adjust  NAve  to  be  an  Odd  Number 
if  NTest==0  %  Even  NAve  Input 

NAve=NAve+l;  %  Adjust  NAve  to  be  Odd 

end; 

for  j=l;Nrow 

%TmpPsd=SigPsd ( j , : ) ; 

SigPsd (j ,: ) =smooth (SigPsd ( j , : ) , NAve) ; 

end 

end 

PSDmax=max (SigPsd. ' ) ; 

PSDf loor=PSDmax/ 10 00 0 ;  %  Set  PSD  Floor  Level  to  -40  dB  Below  Max  Value 

for  j=l:Nrow  %  Threshold  /  Set  PSD  Floor  Values 
for  k=l:Ncol 

if  SigPsd (j , k) <PSDfloor (j ) 

SigPsd ( j , k) =PSDfloor ( j ) ; 

end 

end 

end 

PSDmax=max (SigPsd. ') ;  %  Update  PSDmax 
if  SNorm==l  %  Normalize  PSD  if  Selected 

SigPsd=SigPsd. /repmat ( PSDmax . ' , 1 , Ncol ) ; 

end 

if  NcolIn==l  %  Return  Output  to  Column  Vector 
SigPsd=SigPsd. ' ; 

end 

delf =Fsamp/ (Ncol-1 ) ;  %  Freq  Plot  Step  Size 
Fscale=0 : delf : (Ncol-1) *delf ;  %  Output  PSD  Freq  Scale 

if  SPlot==l  %  Produce  Output  Plots 

disp  (  '  ' ) 

disp  (' Generating  Plot') 
disp  (  '  ' ) 
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if  (Nrow==l | | Ncol==l  ) 

SigPlot= (abs (Sigln) ) .  A2 ; 

SigPlot=SigPlot/max (SigPlot) ; 

PsdPlot=SigPsd; 

else  %  Use  1st  Row  of  Sigln  &  SigPsd  Matrix  for  Plot  Default 
SigPlot= (abs (Sigln (1, : ) ) ) .  A2; 

PsdPlot=SigPsd ( 1,  : )  ; 

end 

figure 

subplot (2 , 1 , 1 )  %  Top  Time  Plot 
tscale=linspace (0,Ncol/Fsamp,Ncol) ; 

SigPlot_dB=10*logl0 (SigPlot) ; 

if  SigVal==l  %  Sigln  is  Real 
plot  (tscale, SigPlot_dB) ; 
grid; 

set  (gca,  'XLim' ,  [0  t_lim] ) ; 
set  (gca,  ' YLim ' ,  [dB_lim  0]  )  ; 
xlabel('Time  [ s ] ' ) 
ylabel ( 'Magnitude ' ) 

else  %  Sigln  is  Complex  . . .  Plot  Amplitude  as  Abs (Sigln) 
plot  (tscale, SigPlot_dB) ; 
grid; 

set  (gca,  ' XLim [0  t_lim] ) ; 
set (gca, ' YLim ', [dB_lim  0] ) ; 
xlabel('Time  [ s ] ' ) 
ylabel ( 'Magnitude ' ) 

end 

save ( ' Time_response ' , ' SigPlot_dB ' , ' tscale ' ) ; 
subplot  (2 , 1 , 2 )  %  Select  Middle  PSD  Figure 

if  SigVal==l  %  Sigln  Real;  Plot  FFTshift{ PSD}  for  -Fs/2  <  f  <  Fs/2 
f pos— ( 0 ; del f ; del f * (ceil (Ncol/2) -1) ) ; 

fneg=  -fliplr(  (delf : delf : delf* (Ncol-ceil (Ncol/2) ) )  ); 

FrqPlot= [ fneg  fpos] ; 

FrqPlot_dB=10*logl0 (fftshift (PsdPlot) ) ; 
plot (FrqPlot, FrqPlot_dB) ;  %  PSD  Shifted 

else  %  Sigln  Complex;  Plot  PSD  for  0  <  f  <  Fsamp 
FrqPlot=0 : delf :  (Ncol-1) *delf; 

FrqPlot_dB=10*logl0 (PsdPlot) ; 

plot (FrqPlot,  FrqPlot_dB) ;  %  FFT  Shifted 


end 

save ( ' PSD_signal ' , ' FrqPlot_dB ' , ' FrqPlot ' ) ; 

%  Adjust  X  &  Y  Axis  Limits 
set  (gca,  'YLim' ,  [dB_lim  0] ) ; 
set  (gca,  'XLim' ,  [0  f_lim] ) ; 
xlabel (' Frequency  [Hz]') 
ylabel ( ' I  FFT  | A  2 ' ) 
grid 

end 

%  ======================  End  SigPSD_V3  Function^ 
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A.2  QMFB  Process  Code 


%  QMFB  PROCESS 
%  Original  file:  qmfb.m 
%  Created  by  :  Phillip  E.  Pace 

%  Modified  by  :  Felipe  Garrido,  CAP.  Chilean  Air  Force. 

%  Last  update  :  02/02/2012. 

%  This  code  is  used  to  create  the  layers  and  the  neccesary  output  data 
%  to  processing  them  according  to  the  desired  time-frequency  plot 
%  This  code  call  the  function  tsinc.m  that  creates  de  coefficients 
%  according  to  the  especified  filter  created  for  Phillip  E.  Pace. 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

clear  all 
close  all 
tic 

%  Files  Loading 

load  preamble_res_10000_25 .mat 
%  --  Input  data  -- 

c0= [preamble_res_10000_25 ] 1 ;  %  Input  signal  padded  with  zeros 

[dl,  cl]  =  tsinc(cO);  %  Sine  Modified  Filter  function 
%  --  Number  of  Layers  to  Compute  -- 
f=c0;  %  Input  signal 

n  =  floor  (log2  (length  (f) )) ;  %  Determine  the  amount  of  layers  from  the  length 

of  the  signal 

N  =  n;  %  number  of  layers 

%  --  Formatting  the  Signal  --  (to  pass  the  signal  through  the  filter  bank) 

I  (1 :2An, 1)  =  f (1 : 2An)  ; 
out  =  I ; 

%  --  Generating  the  Output  Layers  -- 

%  Decompose  the  function 
W  =  waitbar ( 0 Computing  Layers...'); 
for  lay  =  1:N  %  layer 

%disp(lay)  %  Show  what  layer  is  been  generated 

waitbar  (lay/ (N+l) ,W) 

flag  =1;  %  Flag  used  to  set  up  the  columns  in  the  output  matrix 

%  Reshape  the  output  matrix 

[r,  c]  =  size  (out)  ; 

out  =  zeros  (r./2,  c.*2); 

for  i  =  1:2  A (lay-1)  %  column  of  I  (low  to  high) 

[G,H]  =  tsinc  (I  ( : , i) ) ;  %  Evaluate  the  filter  (sine  modified) 

%  over  the  signal  going  by  columns 

%  Setting  the  output  matrix  (Layer)  by  High  Pass  (G)  and  Low  Pass  (H) 
filtering  output 
if  flag 

out ( : , 
out ( : , 

else 

out ( : , 
out ( : , 

end; 

flag  = 

end; 

%  Output  Matrix  and  data  to  work  with 


i .  *2  —  1 )  =  H;  %  Low  Pass  Branch  of  QMFB  Tree 
i.*2)  =  G;  %  High  Pass  Branch  of  QMFB  Tree 

i.*2-l)  =  G;  %  High  Pass  Branch  of  QMFB  Tree 

i . * 2 )  =  H ;  %  Low  Pass  Branch  of  QMFB  Tree 


“f lag; 


%  Change  in  the  flag  value 
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I  =  out;  %  Output  matrix  (signal  already  filtered) 

R  =  real  (I);  %  Real  part  of  the  filtered  signal 

Q  =  imag(I);  %  Imaginary  part  of  the  filtered  signal 

%  Saving  the  data  to  the  same  directory 
save ( [ 1  layer' , int2str (lay) ] ,  'R' ,  '□' ) ; 
end; 

close  (W) 
toe 

%======================End  of  QMFB  Process============================== 

function  [dl,  cl]  =  tsinc  (cO)  %  Sine  Modified  Filter  function 
%  [dl,  cl]  =  tsinc(cO) 

%  Sine  Modified  Filter  function.  It  use  truncated  sine  modified  filter 
%  coefficients  to  decompose  the  column  vector  cO  into  a  (high  frequency) 
%  column  vector  dl  and  (low  frequency)  column  vector  cl. 


%  --  Sine  Modified  Filter  Formula  Application  -- 
N  =  512;  %  Number  of  filters  coefficients  pairs. 

C  =  1.99375872328059;  %  Value  to  get  cross  corr  of  less  than  0.001 
S  =  1.00618488680080;  %  Value  to  get  cross  corr  of  less  than  0.001 
x  =  ( -floor  (N . /2 ):( floor (N . /2 ) -1 ))' ;  %  Vector  coefficients 
h  =  sqrt ( S . /2 )  . *sinc (  (x  +  ,5)./C);  %  Coefficients  filter  formula 
w  =  hamming (N);  %  Hamming  Window 
h  =  w.*h;  %  Low  Pass  Filter  Coefficients 
g  =  h;  %  High  Pass  Filter  Coefficients 

g(2:2:N)  =  -h(2:2:N);  %  Alternating  the  sign  of  the  coefficients 


%  --  Setting  the  Signal  -- 

N  =  length(cO);  %  Length  of  the  signal 

%  (from  qmfb  function  "cO"  is  a  column) 

pad3  =  length  (h) . /2  +  1;  %  pad  with  zeros  to  clear  out  filter 

cO  =  [cO;  zeros (pad3, 1 )] ;  %  padding  with  zeros 


%  Decompose  the  column  vector  cO 
i  =  pad3 : 2 : (N+pad3-2 ) ; 

%  cl  low  frequency  column  vector 
cl  =  filter  (fliplr (h) , 1, cO) ; 
cl  =  cl  (i)  ; 

%  dl  high  frequency  column  vector  dl 
dl  =  filter (fliplr (g) , 1, cO)  ; 
dl  =  dl  (i)  ; 


%  i  will  decimate  by  2 


%  compute  cl 
%  decimate 


%  compute  dl 
%  decimate 
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A. 3  Qn  Matrix  Generation  Code 


Qn  Matrix  generation  and  Plotting 


Original  file 
Created  by 
modified  by 
Last  update 


qmfb .m 

Phillip  E.  Pace 

Capt  Felipe  E.  Garrido  01/04/2012 
02/16/2012 


>ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooc 


%  This  code  allows  to  choose  which  layers  do  you  want  to  plot 
%  and  save  the  matrix  with  useful  data  as  Q  xx.mat,xx=  layer  number 


%  Input  variables  are  the  following: 

%  L=  lower  layer  to  plot  (Min  layer=l) 

%  H=  higher  layer  to  plot  (Max  layer=  N) 

%  fs=  Sampling  frequency  of  the  input  data  (Fsamp_XX) 
%  freq_limit=  max  frequency  in  window 
%  Time  limit=  max  time  in  window 


ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 


clear  all 
close  all 
%%  Warning: 

%  Check  you  are  loading  the  correct  data  file 
load  preamble_res_l 00 00_2 5 . mat ; 


%%  Set  the  appropiate  parameters 


f s=Fsamp_10000; 

L=2  ; 

H=2  4  ; 

f req_limit=9e6 ; 
Time  limit=2 . 5e-5 ; 


%  Sample  frequency  according  to  the 
%  resampled  data 
%  Lower  layer  to  plot 
%  Higher  layer  to  plot 
%  Y  axis  limit 
%  X  axis  limit 


%%  Plotting 
for  yy=L:H 


load  ( [ ' layer ' , num2str (yy) , ' .mat ' ] ) ; 
%  Formating  the  Layer 
layer  =  abs (R  +  li*Q) ;  % 

layer  =  ( layer/max (max ( layer ))). A2 ; % 
M  =  layer;  % 

[m, n] =size  (M)  ;  % 

T=l/fs;  % 

f=linspace ( 0, f s/2 ,  n)  ;  % 

NUMJ3AMPLES  =  m*n;  % 

time=linspace (0,NUM_SAMPLES*T,m) ;  % 

dt= (NUMJ3AMPLES*T) / (m-1) ;  % 

df=  (f s/ (2* (n-1) ) ) ;  % 


absolute  value  of  the  layer 

Normalizing  the  layer 

Set  values  to  dB 

Set  axis  sizes 

Sampling  period 

Setting  the  frequency  axis 

Length  of  the  new  signal 

Setting  the  time  axis 

delta  t 

delta  f 


t_max=ceil (Time_limit/dt) +1 ; 
f_max=ceil ( freq_limit/df ) +1 ; 


new_Q= (M ( 1 : ( tmax ) , 1 : ( f_max ))).'; 

new_layer= (R ( 1 : ( t_max) , 1 : ( f_max) )  +  li*Q ( 1 : ( t_max) , 1 : ( f_max) ) ) . ' 
n  t=time ( 1 : t_max) ; 
n_f=f (1 : f_max) ; 

Q  dB=10*logl0 (new  Q.A2); 
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o\°  o\°  o\°  o\°  o\°  o\° 


clear  layer; 


%  clear  current  layer  matrix  from  the  workspace 


%  Color  Plots 
figure; 

%pcolor (time ( 1 : tmax) , f ( 1 : f_max) , Q_dB ' ) 
pcolor (n_t, n_f , Q_dB) 

title ( [ ' "Time  -  Frequency  Plane  Layer  ' , . . . 

num2str (yy) , '  with  \Deltat= ' , num2str (dt) , . . . 

'  and  \Deltaf=' , num2str (df ) 
xlabel('Time  [s] '), ylabel (' Frequency  [Hz]') 
colorbar, colormap (jet) 
caxis ( [-30, 0] ) 
shading  interp 

figure; 

%pcolor (time ( 1 : t_max) , f ( 1 : f_max) , Q_dB ' ) 
pcolor (n_t, n_f , Q_dB) 

title ( [ ' "Time  -  Frequency  Plane  Layer  ' , . . . 

num2str (yy) , 1  with  \Deltat= ' , num2str (dt) , . . . 

'  and  \Deltaf= num2str (df ),'"'] ) 
xlabel('Time  [sec] '), ylabel (' Frequency  [Hz]') 
colorbar, colormap (jet) 
caxis ( [-30, 0] ) 
shading  flat 

save ( [ ' Q _ ' , int2str (yy) ] , ' Q_dB ' , ' new_Q ' , ' time ' , ' f ' , ' dt ' , ' df ' , . . . 

' f s ' , ' tmax ' , ' f_max ' , ' new_layer ' ) ; 

end 
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A.4  Input/Output  Frequency-Time  (F-T)  Analysis  Code 


%  TIME  AND  FREQUENCY  ANALYSIS 

%  Created  by  :  Felipe  E.  Garrido,  CAP.  Chilean  Air  Force 
%  Last  update  :  02/16/2012 

"o 

%  This  code  computes  and  save  Normalized  Average  Time 
%  and  Normalized  Average  frequency  vector  for  a  given  layer 
"6 

'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O 

%%  Warning 

%  1.-  The  file  Layer_plotting,  must  be  run  before  running  this  code  to 
%  generate  a  given  Q_XX.mat  file. 

'o 

oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 
"6  "6 

%  Input  variables  are  the  following: 

%  layer  =  number  of  the  layer  you  want  to  plot  (Min  layer=l) 

%  sel  =  variable  to  choose  between  time  or  frequency 
%  Choose: 

%  Time  plot  ===>  sel=l 

%  Frequency  plot  ===>  sel=2 
"6 

oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

clear  all 
close  all 

layer=13;  %  Choose  the  layer  to  plot 

load  (['□_’ ,num2str (layer) , ' .mat' ] ) ; 


sel=l;  %  Time  plot=l.  Frequency  plot=2 

%  n  t=time(l:t  max) ; 

%  n_f =f ( 1 : f_max) ; 

%  Q_dB=10*logl0 (new_Q . A2 ) ; 

if  sel==l; 

t=mean (new_Q, 1) ; 
t2=t/max (t) ; 

t2_dB=10*logl0 (t2 (1 : t  max) ) ; 

tt=time (1 :t_max) ; 

figure; 

plot (tt, t2_dB) 

%  title ( [ ' "Normalized  Average  Time  plot  for  Layer  ' , . . . 

%  num2str (layer) , '  with  \Deltat= ' , num2str (dt) , ' " ' ] ) 

%xlabel ( ' Time  [ s ] ' ) 
ylabel ( 'dB' ) 
ylim ( [-30  0] ) 
xlim([0  Time_limit] ) 
grid  on 

save ( [ ' ave_t_ ' , int2str (layer ) ] , ' t2_dB ' , ' tt ' ) ; 

elseif  sel==2; 

nf=mean (new_Q, 2 ) ; 
f 2=nf /max (nf ) ; 
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f2_dB=10*logl0 ( f 2 (1 : f_max) ) ; 

w=f (1 : f_max) ; 

figure; 

plot ( f2_dB, w) 

%  title ( [ ' "Normalized  Average  frequency  plot  for  Layer  ' , . . . 

%  num2str (layer) , '  with  \Deltaf= num2str (df ) 

%ylabel (' Frequency  [Hz]') 
xlabel ( 'dB' ) 
xlim ( [-30  0] ) 
ylim([0  freq_limit] ) 
grid  on 

save ( [ ' ave_f_ ' , int2str (layer ) ] , ' f  2_dB '  , '  w '  )  ; 

else 

disp ( 'ERROR: ' ) 

disp (' Choose :  l=Plot  Average  Time  or  2=Plot  Average  Frequency') 

end 

%=====================End  of  Time/Frequency  analysis  code========== 


%  OVERLAY  FREQUENCY  RESPONSE  PLOT 
"6 

%  Created  by  :  Felipe  E.  Garrido,  CAP.  Chilean  Air  Force 
%  Last  update  :  02/16/2012 

%  This  code  computes  overlay  frequency  response  according  to  input 
%  signal  and  two  chosen  layers 


oooooooooc 


lOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 


%  Input  variables  are  the  following: 

%  PSD_signal . mat  =  initial  signal  frequency  response 
%  ave_f_20.mat  =  single  burst  frequency  response 
%  ave_f_19.mat  =  burst  integration  frequency  response 

'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O 


clear  all 
close  all 

load  PSD_signal .mat 
figure 

plot  (FrqPlot, FrqPlot_dB, 'b ' ) 
hold  on 

load  ave_f_20.mat 
plot (w, f2_dB, ' r ' ) 

hold  on 

load  ave_f_19.mat 
plot (w, f2_dB, 'k') 

grid  on 

legend ( 'Input' , 'Single  Burst' , 'Burst  Integration' , 'location' , 'South' ) 

xlim( [0  2e6] ) 

ylim ( [-30  0]  ) 

xlabel (' Frequency  [Hz]') 

ylabel ( 'Magnitude  [dB] ' ) 

%=====================End  of  Overlay  FREQUENCY  response============== 
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%  OVERLAY  TIME  RESPONSE  PLOT 


%  Created  by  :  Felipe  E.  Garrido,  CAP.  Chilean  Air  Force 
%  Last  update  :  02/16/2012 
"6 

%  This  code  computes  overlay  frequency  response  according  to  input 
%  signal  and  two  chosen  layers 


ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

%  Input  variables  are  the  following: 

%  Time  response. mat  =  initial  signal  frequency  response 
%  ave_t_13.mat  =  single  burst  frequency  response 
%  ave  t  12. mat  =  burst  integration  frequency  response 


'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O 


clear  all 
close  all 

load  Time_response . mat 
figure 

plot  (tscale, SigPlot_dB, 'b ' ) 
hold  on 

load  ave_t_13.mat 
plot (tt, t2_dB, 'r') 

hold  on 

load  ave_t_12.mat 
plot (tt, t2_dB, ' k' ) 

grid  on 

legend ( ' Input ' ,  ' Single  Burst ' ,  ' Burst  Integration  1 ,  1  location  1 ,  ' south ' ) 

xlim ( [0  0.1e-3] ) 

ylim ( [-30  0] ) 

xlabel('Time  [ s ]  ' ) 

ylabel ( 'Magnitude  [dB] ' ) 

%=====================End  of  Overlay  TIME  response=================== 
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Frequency  [Hz] 


APPENDIX  B  -  QMFB  Outputs  for  802.1  la  Wi-Fi  Preamble  Signal 


This  appendix  presents  a  complete  set  of  available  layer  outputs  (Qn)  for  the  802.1  la 
signal  analyzed  in  Section  4.2  according  to  computed  QMFB  outputs. 
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QMFB  Layer  Q? 


QMFB  Layer  Q3 
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Frequency  [Hz]  Frequency  [Hz]  Frequency  [Hz] 


QMFB  Layer  Q(, 


QMFB  Layer  Qj 


Time  [s]  x  ^g'5  Time  [s] 


QMFB  Layer  Q$ 


QMFB  Layer  Qg 


QMFB  Layer  Qm 


QMFB  Layer  Qu 
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QMFB  Layer  Qn 


QMFB  Layer  Qu 


QMFB  Layer  Qm 


e 


QMFB  Layer  Qis 


QMFB  Layer  QI6 


6 


Time  [s] 


QMFB  Layer  Qn 
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Frequency  [Hz]  Frequency  [Hz]  Frequency  [Hz] 


6 


QMFB  Layer  Q19 


QMFB  Layer  Q20 


QMFB  Layer  Q21 


QMFB  Layer  Q22 


QMFB  Layer  Q23 
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6 


QMFB  Layer  Q24 
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APPENDIX  C.  Time  (At)  and  Frequency  (Af)  Resolution  Tables 


C.l  802.1  la  Wi-Fi  Signal:  Preamble  Only  (Section  4.2.1) 


Sampling 
Frequency 
fs  [Hz] 

fs  =  5.94xl010 

fs=  1.19X1011 

fs  =  2.38x10“ 

Layer 

Af 

At 

Af 

At 

Af 

At 

# 

[Hz] 

[s] 

[Hz] 

[s] 

[Hz] 

[s] 

6 

4.71xl08 

1.08xl0'9 

9.43x10s 

5.39xlO'10 

1.89xl09 

2.70xl0'10 

7 

2.33x10s 

2. 16x1  O'9 

4.68x10s 

1.08x1  O'9 

9.35xl08 

5.39xlO'10 

8 

1.16xl08 

4.31xl0'9 

2.33x10s 

2.16xl0'9 

4.66xl08 

1.08xl0‘9 

9 

5.81xl07 

8.62xl0'9 

1.16x10s 

4.31xl0'9 

2.32xl08 

2.16xl0‘9 

10 

2.90xl07 

1.73xl0'8 

5.80xl07 

8.62x1  O'9 

1.16xl08 

4.31xl0‘9 

11 

1.45xl07 

3. 45x1  O'8 

2.90xl07 

1.73xl0'8 

5.80xl07 

8.62x1 0‘9 

12 

7.25xl06 

6.90xl0'8 

1.45xl07 

3.45xl0'8 

2.90xl07 

1.73xl0‘8 

13 

3.62xl06 

1.38xl0'7 

7.25xl06 

6.90xl0'8 

1.45xl07 

3.45xl0'8 

14 

1.81xl06 

2.76xl0'7 

3.62xl06 

1.38xl0'7 

7.25xl06 

6.90xl0'8 

15 

9.06xl05 

5.52xl0'7 

1.81xl06 

2.76xl0'7 

3.62xl06 

1.38xlO'7 

16 

4.53x10s 

1 . 1  lxl  0"6 

9.06xl05 

5.53xl0'7 

1.81xl06 

2.77xl0'7 

17 

2.27xl05 

2.22x1  O'6 

4.53x10s 

1 . 1  lxl  0‘6 

9.06xl05 

5.56xl0'7 

18 

1.13xl05 

4.45xl0'6 

2.26xl05 

2.23xl0'6 

4.53x10s 

1 . 1  lxl  0‘6 

19 

5.66xl04 

8.97xl0'6 

1.13xl05 

4.49x1  O'6 

2.27xl05 

2.25xl0‘6 

20 

2.83xl04 

1.82x1  O'5 

5.66xl04 

9.12xl0'6 

1.13xl05 

4.56xl0‘6 

21 

1.42xl04 

3. 77x1  O'5 

2.83xl04 

1.88xl0'5 

5.66xl04 

9.42x1  O’6 
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C.2  802.1  la  Wi-Fi  Signal:  Complete  Burst  (Section  4.2.2) 


Sampling 

Frequency 

fs  [Hz] 

f  =  2.38xl010 

fs  =  5.94  xlO10 

fs=  1.31  xlO11 

Layer 

Af 

At 

Af 

At 

Af 

At 

# 

[Hz] 

[s] 

[Hz] 

[s] 

[Hz] 

[s] 

6 

1.88x10s 

2.69x1  O'9 

4.71x10s 

1.08x1  O'9 

1.04xl09 

4.90xl0'10 

7 

9.35xl07 

5.39xl0'9 

2.33x10s 

2. 16x1  O'9 

5.14x10s 

9.80xl0'10 

8 

4.66xl07 

1.08xl0's 

1.16x10s 

4.31xl0'9 

2.56x10s 

1.96x1  O'9 

9 

2.32xl07 

2.16xl0's 

5.81xl07 

8.62xl0'9 

1.27x10s 

3. 92x1  O'9 

10 

1.15xl07 

4.31xl0's 

2.90xl07 

1.73xl0'8 

6.38xl07 

7.84xl0'9 

11 

5.80xl06 

8.63xl0's 

1.45xl07 

3. 45x1  O'8 

3.19xl07 

1.57xl0's 

12 

2.89xl06 

1.73xl0'7 

7.25xl06 

6.90xl0'8 

1.59xl07 

3.14xl0's 

13 

1.44xl06 

3. 45x1  O'7 

3.62xl06 

1.38xlO'7 

7.97xl06 

6.27xl0's 

14 

7.25x10s 

6.93xl0'7 

1.81xl06 

2.76xl0'7 

3.99xl06 

1.25xl0'7 

15 

3.62xl06 

1.39xl0'6 

9.06xl05 

5.52xl0'7 

1.99xl06 

2.51xl0'7 

16 

1.81xl05 

2.80xl0'6 

4.53x10s 

1.1  lxl  0"6 

9.97xl06 

5.03xl0'7 

17 

9.06xl04 

5. 70x1  O'6 

2.27xl05 

2.22x1  O'6 

4.98xl05 

l.OlxlO'6 

18 

4.53xl04 

1.18xl0'5 

1.13xl05 

4.45xl0'6 

2.49xl05 

2.02xl0'6 

19 

2.26xl04 

2.52xl0'5 

5.66xl04 

8.97xl0'6 

1.25x10s 

4.08xl0'6 

20 

1.13xl04 

5.89xl0'5 

2.83xl04 

1.82x1  O'5 

6.23xl04 

8.29xl0'6 

21 

5.63xl03 

1.21xl0'4 

1.42xl04 

3. 77x1  O'5 

3.1 lxlO4 

1.71xl0'5 
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C.3  802. 16e  WiMAX  Signal 


Sampling 

Frequency 

fs  [Hz] 

Both  Modes 
fs  =  1.19xl09 

Data-Only 
fs=  1.18x10*° 

Range-Only 
fs  =  3.56  xlO*° 

Layer 

Af 

At 

Af 

At 

Af 

At 

# 

[Hz] 

[s] 

[Hz] 

[s] 

[Hz] 

[s] 

6 

9.42xl06 

5.39xl0'8 

9.42xl07 

5. 39x1  O'9 

2.83xl08 

1.80xl0'9 

7 

4.68xl06 

1.08xl0'7 

4.68xl07 

1.08xl0'8 

1.40xl08 

3.59xl0'9 

8 

2.32xl06 

2. 16x1  O'7 

3.33xl07 

2. 16x1  O'8 

6.99xl07 

7. 19x1  O'9 

9 

1.15xl06 

4.31xl0'7 

1.16xl07 

4.31xl0'8 

3.49xl07 

1.44x1  O'8 

10 

5.80x10s 

8.62xl0'7 

5.80xl06 

8.62xl0'8 

1.74xl07 

2.87xl0'8 

11 

2.90xl05 

1.72x1  O'6 

2.90xl06 

1.72x1  O'7 

8.70xl06 

5.75xl0'8 

12 

1.45x10s 

3.45xl0'6 

1.45xl06 

3.45xl0'7 

4.34xl06 

1.15xlO'7 

13 

7.24xl04 

6.90xl0'6 

7.25x10s 

6.90xl0'7 

2.16xl06 

2.31xl0'7 

14 

3.62xl04 

1.38xl0'5 

3.62xl05 

1.38xl0'6 

1.09xl06 

4.60xl0'7 

15 

1.81xl04 

2.76xl0'5 

1.81xl05 

2.76xl0'6 

5.44x10s 

9.21xl0'7 

16 

9.06xl03 

5.53xl0'5 

9.06xl04 

5.53xl0'6 

2.71xl05 

1.84x1  O'6 

17 

4.53xl03 

1 . 1  lxl  0"4 

4.53xl04 

l.llxlO'5 

1.36xl05 

3. 69x1  O'6 

18 

2.27xl03 

2.33xl0'4 

2.27xl04 

2.22x1  O'5 

6.80xl04 

7.42x1  O'6 

19 

1.13xl02 

4.75xl0'4 

1.13xl04 

4.45xl0'5 

3.39xl04 

1.50xl0'5 

20 

5.63x10* 

1.02x1  O'3 

5.66xl03 

8.97xl0'5 

1.69xl04 

3.03xl0'5 

21 

2.84x10* 

2.1  lxlO'3 

2.83xl03 

1.82x1  O'4 

8.58xl03 

6.04xl0'5 
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